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A  general  description  Is  given  of  the  computer  program  which  simulates 
the  performance  of  an  anti-radiation  homing  weapon  talcing  into  account  the 
properties  of  the  seeker,  assumed  to  be  body-referenced,  and  the  air  vehicle. 
For  the  purpose  of  Illustration,  sample  calculations  have  assumed  a  mlnldrone 
as  the  air  vehicle,  although  the  simulation  has  a  wide  range  of  applications 
to  a  variety  of  platforms.  The  program  can  simulate  terminal  homing  in  the 
situations  where  multiple,  co-channel  CW  emitters  are  deployed.  The  three 
constituent  simulation  blocks  (signal  in  space,  seeker,  and  tactical  vehicle) 
are  discussed  in  detail,  as  well  as  models  for  sensor  errors  and  wind  effects. 
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1.0  INTRODUCTION  AND  OVERVIEW 

A  coaputer  program  has  been  developed  for  Che  purpose  of  evaluation  of  an 
emitter  homing  seeker  developed  by  Group  44.  The  details  of  the  computations 
required  to  model  the  electromagnetic  environment,  and  seeker  and  vehicle 
characteristics  are  given  in  various  technical  reports.  In  particular,  the 
model  for  the  electromagnetic  environment  is  described  in  [1].  The  simulation 
can  be  considered  to  consist  of  four  parts: 

(1)  An  initialization  section 

(2)  A  signal  in  space  section,  modelling  target  attributes  and 
various  multipath  effects  such  as  specular  and  diffuse  re¬ 
flections  off  the  ground,  obstacles  and  shadowing 

(3)  A  seeker  module,  modelling  antenna  pattern,  signal  proces¬ 
sing  algorithm,  and  tracker  and  guidance  logic 

(4)  A  control  and  aerodynamics  block,  modelling  vehicle  autopi¬ 
lot  and  aerodynamics. 

The  driver  for  these  parts  is  Program  SJCMLST.  For  illustrative  pur¬ 
poses,  the  vehicle  and  seeker  that  are  considered  are  the  XBQM-106  minidrone 
and  the  enhanced  interferometer. 

A  generalized  block  diagram  is  given  in  Fig.  1.1.  The  overall  simulation 
program  sequence  is  illustrated  in  Fig.  1.2.  This  simulation  sequence  diagram 
is  typical  of  the  diagrams  which  will  be  used  to  describe  the  structure  of  the 
program.  The  blocks  in  the  diagram  represent  programs,  subprograms  or  func¬ 
tions,  and  the  arrows  Indicate  the  direction  of  the  call  and  return.  As  may 
also  be  seen  on  this  diagram,  an  indentation  in  the  blocks  is  used  to  empha¬ 
size  the  heirachy  of  calling  dependence. 

The  simulation  components  are  covered  in  Section  2.0  through  5.0.  Sec¬ 
tion  6.0  discusses  the  EXEC  files  which  are  used  to  run  the  program  on  the 
Amdahl  480  computer  at  Lincoln  Laboratory.  Section  7.0  presents  an  example  of 
a  typical  simulation  output. 

Appendices  A  through  D  give  further  details  of  the  individual  blocks  il¬ 
lustrated  in  the  program  sequence  diagrams,  including  input  and  output  as  well 
as  definitions  of  the  variables  used.  These  appendices  are  most  useful  when 
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Fig.  1.1.  Generalized  simulation  block  diagram. 
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SJCMLST 


Fig.  1.2.  Overall  simulation  program  processing. 
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considered  In  conjunction  with  the  program  listings  which  accompany  this  re¬ 
ports  The  programs  themselves  are  written  primarily  in  FORTRAN  and  structured 
FORTRAN  (IFTRAN). 


2.0  INITIALIZATION  MODULE 

The  program  sequence  for  the  initialization  module  is  shown  in  Fig.  2.1. 
As  noted  in  Section  1.0,  the  call  to  the  module  is  from  SJCMLST,  the  master 
driver  of  the  simulation.  The  purpose  of  this  module  is  to  both  initialize 
the  appropriate  variables  as  well  as  to  read  in  information  from  the  EXEC, 
which  specifies  the  seeker  angle  processing  options  to  be  considered  (such  as 
the  enhanced  interferometer  or  amplitude  monopulse;.  The  input  information 
specified  by  the  user  in  block  data  file  BDUHFM  and  the  EXEC  is  written  out  to 
the  terminal.  Data  specified  in  BDUHFM  is  manipulated  to  set  up  arrays  which 
are  used  in  the  subroutines  modelling  the  reflections  off  the  ground,  build¬ 
ings,  and  aircraft.  The  driver  SJCMLST  considers  each  of  the  receiver  options 
in  turn  and  determines  the  vehicle  trajectory  and  seeker  performance  for  each 
of  the  direction  finder  algorithms  considered.  A  number  of  criteria  are  em¬ 
ployed  to  determine  when  to  switch  from  consideration  of  one  seeker  to  anoth¬ 
er.  Following  the  initialization  steps,  multipath  attributes,  seeker  calcula¬ 
tions,  and  flight  dynamics  calculations  are  sequenced  as  shown  in  Fig.  1.2. 
This  section  will  consider  in  a  general  way  the  flow  of  the  initialization, 
while  Appendix  A  will  given  additional  details  of  the  subroutines,  along  with 
definitions  of  the  variables  used  in  the  program  listings  which  accompany  this 
report. 

Subroutine  PTIME  is  called  initially  to  provide  a  beginning  time  which  is 
subsequently  used  in  calculations  to  determine  the  CPU  time  needed  for  multi- 
path,  seeker  and  flight  dynamics  procedures.  MISINT  is  included  for  the  ini¬ 
tialization  of  two  logical  variables  associated  with  the  Maverick  missile  and 
is  not  applicable  for  the  current  simulation  of  the  XBQM-106.  The  block  la¬ 
belled  (1)  indicates  that  several  initializations  are  performed  and  emitter 
frequencies  calculated  with  lines  of  code  in  the  driver  itself. 

The  subroutine  READIT,  which  c-'lls  a  number  of  other  subroutines,  reads 
the  number  which  indexes  the  run  and  the  title  which  labels  the  run.  These 
are  specified  in  the  EXEC.  The  EXEC  input  is  also  used  to  determine  whether 
specular  ground  reflections  will  be  modelled  as  well  as  to  provide  initial 
conditions  for  the  vehicle  state  vector  (orientation,  position,  velocity). 
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The  subroutine  RORI  reads  the  number  used  as  a  label  to  Index  the  simula¬ 
tion  run.  The  subroutine  YESNO  receives  the  user's  response  as  to  whether 
ground  reflections  are  to  be  modelled.  This  subroutine  calls  RDRC  which  ac¬ 
tually  reads  the  input  from  the  EXEC.  The  subroutine  PUTCHR  is  used  in  the 
validation  check  of  the  user's  response  to  the  input  questions.  The  subrou¬ 
tine  WRTYN  writes  'yes'  or  'no'  to  the  terminal  regarding  the  modelling  of 
specular  ground  reflections. 

The  subroutine  HOJRDT  reads  which  receivers  are  to  be  considered  for  the 
simulation  run .  Capability  exists  to  model  a  wide  variety  of  candidate 
antenna/receiver/signal  processing  systems  for  direction  finding.  It  also  ac¬ 
cepts  as  input  the  vehicle  state  vector  with  respect  to  the  reference  Inertial 
coordinate  system.  To  accomplish  these  tasks  a  number  of  subroutines  are  in¬ 
voked.  GETRCV ,  RDRCV ,  AND  PUTCHR  perform  the  selection  of  the  receiver  op¬ 
tion.  RDXYZ  reads  the  specified  position  of  the  vehicle  in  space  while 
INTPOS,  which  calls  INVERT,  uses  the  initial  specification  of  velocity  and 
boresight  orientation  with  respect  to  Emitter  1  to  determine  the  vehicle  bore- 
sight  vector  with  respect  to  Inertial  space. 

The  block  labelled  (2)  indicates  that  there  are  lines  of  code  in  the 
driver  itself  which  checks  the  number  of  obstacles  (buildings  or  aircraft), 
emitters  and  blinking  Information  to  insure  that  the  number  of  input  quanti¬ 
ties  remain  within  the  allowed  dimensions. 

TITLPG  writes  out  some  of  the  parameters  associated  with  the  emitters, 
such  as  position,  wave  length,  and  polarization  on  the  first  page  of  the  out¬ 
put.  On  the  second  page  of  the  output,  parameters  associated  with  the  multi- 
path  are  written  out. 

flXMURT  writes  out  the  terminal  blinking  emitter  parameters,  initial  con¬ 
ditions  of  the  vehicle,  and  parameters  associated  with  the  seeker. 

The  subroutine  INIT,  which,  in  turn  calls  a  number  of  other  subroutines, 
initializes  the  common  areas  for  the  signal  in  space  routines  based  on  the  in¬ 
put  data  contained  in  the  block  data  file  BDUHFM. 
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The  subroutine  SURFSW  modifies,  if  necessary,  an  array  containing  the 
corners  of  the  rectangles  and  triangles  describing  the  ground  plane.  Other 
subroutines  that  are  subsequently  called  expect  the  corners  to  be  entered  in  a 
certain  order.  This  subroutine  checks  the  appropriate  array  and,  if  neces¬ 
sary,  modifies  it  to  insure  that  the  correct  order  is  maintained. 

The  subroutines  VSUB,  ORIENT,  DOT,  and  NORML  along  with  its  called  rou¬ 
tines  CROX  and  NORMLZ  are  a  series  of  vector  manipulation  subroutines  needed 
for  INIT.  The  subroutine  ACTYPE  specifies  certain  parameters  for  aircraft. 
For  the  simulations  of  scenarios  which  do  not  contain  reflections  off  aircraft 
this  subroutine  has  no  effect  on  the  results. 

The  call  to  the  subroutine  SCLHOJ  establishes  scaling  parameters  for  the 
plotting  programs. 

The  block  labelled  (3)  represents  lines  of  code  in  the  driver  $JCMLST  for 
several  initializations  after  which  calls  are  made  to  the  other  main  modules 
of  the  simulation  to  determine  the  time  history  of  the  seeker/ vehicle  system 
in  the  tactical  scenario  specified  in  the  block  data  file  BDUHFM. 

Return  to  the  driver  is  cycled  with  calls  to  the  subroutine  FLIGHT  to 
determine  when  a  given  receiver  loop  is  to  be  ended  and  final  miss  distances 
are  to  be  determined.  The  criteria  used  for  this  are  as  follows: 

(a)  the  altitude  of  the  vehicle  is  less  than  0.1  ft. 

(b)  the  altitude  of  the  vehicle  is  less  than  the  height  of  a 
particular  emitter  and  the  distance  from  the  emitter  is  less 
than  50  ft. 

(c)  the  index  which  specifies  the  flight  time  exceeds  a  speci¬ 
fied  value. 

The  subroutines  used  in  FLIGHT  to  accomplish  these  tasks  are: 

XDIV  performs  the  normalization  of  a  three  component  vector 

MINPOS  determines  the  miss  distance  to  emitter  IXMTR  and  calls 

EMCASE  which  determines  the  coordinates  of  a  vector  from  the 
vehicle  to  the  emitter. 
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which  determines  Che  alee  dl stance  using  a  point  of 
closest  epproach  analysis. 


The  subroutine  MOVE  Is  used  to  load  emitter  parameters  Into  a  labelled 
common  area. 

The  block  labelled  (4)  Indicates  lines  In  the  driver  used  to  make  a  num¬ 
ber  of  Initializations  after  which  the  subroutine  IM1TR  is  called  to  compute 
the  receiver  position  and  velocity  coordinates  in  a  reference  frame  of  equiva¬ 
lent  vertical  aircraft  fuselages. 

The  subroutine  XNITXR  computes  the  relative  vector  from  the  transmitter 
to  the  receiver,  the  wave  number  and  the  wave  number  squared.  This  subroutine 
calls  the  vector  subroutines  VSUB,  VNORM  and  VNORM. 

The  final  block  designated  (5)  represents  lines  of  code  in  the  driver  to 
calculate  the  ratio  of  the  distance  from  the  receiver  to  Emitter  1  along  with 
the  associated  time  delay. 
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3.0  SIGNAL  IN  SPACE  MODULE 

The  signal  In  space  module  characterizes  the  electromagnetic  environment 
as  seen  by  the  receive  antenna.  As  shown  in  Fig.  3.1  this  is  done  in  a  multi¬ 
emitter  scenario  by  considering  each  of  the  emitters  in  turn  and  computing  the 
nultlpath  parameters  associated  with  the  direct  wave  from  emitter  to  the  re¬ 
ceiver  in  the  vehicle  along  with  the  specular  reflections  off  the  ground  and 
other  obstructions  such  as  buildings  and  other  aircraft  in  the  vicinity  of  the 
emitters.  It  is  to  be  noted  that  different  bounce  paths  are  assumed  for  re¬ 
flections  off  the  various  obstructions.  Figure  3.2,  which  describes  the  over¬ 
all  program  flow  of  the  module,  shows  that  each  of  the  major  submodules,  or 
subroutines,  are  called  by  Che  master  similation  driver,  $JCMLST. 

Subsection  3.1  discusses  the  sequencing  of  the  subroutines  which  are  used 
in  the  direct  wave  computation,  while  subsections  3.2,  3.3,  and  3.4  review  the 
subroutine  used  for  the  specular  ground,  building,  and  airplane  reflections, 
respectively.  Subsection  3.5  discusses  the  parameters  required  for  multipath 
computations.  Appendix  B  provides  additional  details  of  the  input  and  output, 
as  well  as  the  definition  of  variables  given  in  the  program  listings  which  ac¬ 
company  this  report. 

In  most  of  the  simulation  runs  done  for  the  EHT  studies,  building  and 
aircraft  reflections  as  well  as  shadowing  due  to  ground  humps  were  not  mod¬ 
elled  but  are  included  in  this  description  for  completeness. 

3.1  Direct  Wave  Computation 

The  program  sequence  for  the  direct  wave  computation  is  shown  in  Fig. 
3.3.  This  part  of  the  signal  in  space  program  determines  the  multipath  param¬ 
eters  associated  with  the  direct  wave  from  the  emitters  including  any  shadow¬ 
ing  effects  associated  with  obstructions.  The  first  step  in  the  procedure  is 
to  Invoke  the  subroutine  CHKMLT  which  checks  the  number  of  multipath  compo¬ 
nents,  and,  if  there  are  more  than  a  specified  number,  stops  execution  of  the 
program.  Next,  the  subroutine  DIRWV  computes  the  multipath  parameters  associ¬ 
ated  with  the  direct  wave,  Including: 
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ENTER  MULTIPATH  COMPUTATION  PROCEDURE 


EXIT  MULTIPATH  COMPUTATION  PROCEDURE 

Fig.  3.1.  Multipath  computation  procedure 
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Fig.  3.2  Overview  of  program  processing  for  wave  in  space 
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Fig.  3.3  Prograa  processing 


*  The  Doppler  shift 

*  The  planar  angles  which  specify  the  direction  of  the  inci¬ 

dent  radiation  with  respect  to  the  receive  antenna 

*  The  planar  angles  which  specify  to  the  location  of  the  re¬ 

ceiver  relative  to  the  transmitter 

In  accomplishing  these  tasks ,  several  subroutines  are  called: 

RCVANG  which  calculates  planar  angles  in  the  receiver  refer¬ 
ence  system  given  the  coordinates  relative  to  the  ref¬ 
erence  inertial  coordinates. 

VNORMC  which  coa$utes  the  length  of  a  three-component  vector 
DOTC  which  determines  the  dot  product  of  two  vectors 

ATAN  which  finds  the  arc  tangent 

Next,  the  subroutine  SHDHMJ  models  the  shadowing  effect  due  to  circular 
cylindrical  humps  and  determines  the  multipath  parameters  associated  with 
these  effects.  Figure  3.4  shows  the  standard  coordinate  system  used  in  the 
computations  done  in  this  subroutine.  As  noted  in  Fig.  3.3,  subroutine  SHDHMJ 
calls  a  number  of  other  subroutines,  the  major  two  being  CIRCSH  and  CSPK. 
Subroutine  CIRCSH  confutes  the  parameters  needed  to  determine  the  shadowing 
effect  due  to  a  circular  cylindrical  hump.  CIRCSH  calls,  in  turn,  GZ  which 
computes  the  complex  function  required  to  determine  the  shadowing  of  the  hump 
and  HMPROT  which  transforms  inertially  reference  coordinates  of  a  vector  into 
a  coordinate  system  based  on  the  vector  from  the  transmitter  to  the  receiver. 
The  subroutine  CSPR  determines  the  Fresnel  Integral. 

The  subroutine  SHDHMJ  will  .lot  be  called  if  the  parameter  IHUMP  is  set  to 
zero  (its  default  value)  in  the  block  data  file  BDUHFM. 

3.2  Ground  Reflection  Computation 

Figure  3.5  outlines  the  sequence  for  ground  reflection  computation.  As 
for  the  direct  wave  calculation,  the  subroutine  CHKMLT  is  first  called  to  in¬ 
sure  that  the  number  of  multipath  components  does  not  exceed  a  specified 
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Fig.  3.5.  Program  processing  for  ground  reflection 
computation 
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value.  If  specular  ground  reflections  are  modelled  (as  in  most  simulation 

v 

cases),  subroutine  GREFC  is  then  called.  If  not,  then  the  parameters  associ¬ 
ated  with  multipath  are  set  to  zero.  If  the  parameter  ISPGRD  is  specified  in 
block  data  BDUHFM  as  l  (again,  its  default  value)  a  simplified  method  is  used 
to  find  the  oultipath  parameters  in  subroutine  FRERFS.  On  the  other  hand,  if 
ISPGRD  is  zero,  the  Fresnel  integral  is  calculated  using  a  two-dimensional  ap¬ 
proximation.  In  the  latter  case,  the  ground  plane  can  be  represented  by  a  set 
of  arbitrarily  oriented  rectangles  and  triangles,  each  with  its  own  character¬ 
istic  surface  roughness  coefficient  and  dielectric  constant.  A  flow  chart 
describing  the  computation  of  multipath  parameters  is  given  in  Fig.  3.6.  As 
may  be  noted  in  this  diagram  specular  reflection  parameters  are  computed  for 
each  region  or  element.  The  simplified  method  bypass  (for  ISPGRD  «  1)  which 
uses  default  values  of  the  surface  roughness  and  dielectric  constant  may  also 
be  seen  in  the  diagram. 

Multipath  parameters  are  adjusted  by  modelling  the  transmit  antenna  pat¬ 
tern  in  subroutines  WGTAMP  and  GXMTR.  The  multipath  parameters  are  loaded  in¬ 
to  labelled  common  areas  MULTI  and  MULTIP  by  means  of  subroutine  SETVAL.  Sub¬ 
routine  DIGLOD  transfers  the  parameters  associated  with  the  maximum  amplitude 
multipath  component  into  labelled  common  area  DIAGP  after  being  modified  by 
the  receiver  antenna  pattern  which  is  characterized  in  GAIN. 

Appendix  B  supplies  additional  information  on  the  subroutine  discussed 
and  shown  in  Fig.  3.5. 

3.3  Building  Reflection  Computations 

The  sequence  of  computation  for  building  reflections  is  shown  in  Fig. 
3.7.  When  buildings  are  included  in  a  scenario,  calculations  are  made  for 
each  building  specified  considering  four  different  bounce  paths  as  follows: 
transraitter-ground-building-receiver 
transmitter-building-ground-receiver 
transmi tter-ground-building-ground-receiver 
transmitter-building-receiver 
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Fig.  3.6.  Flow  chart  for  computation  of  multipath 
parameters  for  specular  ground  reflection. 
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As  for  Che  other  computation  segments  In  this  module,  CHKMLT  la  called  flrBt 
to  Insure  that  Che  number  of  multipaths  does  not  exceed  the  specified  number. 
Next,  the  subroutine  BREFC  is  called  to  calculate  the  parameters  associated 
with  reflection  off  buildings  such  as  amplitude,  phase,  time  delay,  and 
Doppler  shift.  BREFC,  in  turn,  calls  a  series  of  other  subroutines  as  shown 
in  Fig.  3.7.  Subroutines  BCOMP,  BCOORD,  BCOORT  are  coordinate  transformation 
subroutines  and  VSUB ,  VNORM,  VNORMC  are  vector  operation  subroutines.  RCVANG 
is  a  subroutine  that  is  used  to  calculate  planar  angles  relative  to  vehicle 
axes  given  an  inertially  referenced  vector  from  the  vehicle  to  a  point  in 
space.  The  subroutine  CSPR,  discussed  in  Section  3.1,  determines  the  Fresnel 
integral,  while  POLZN  computes  the  equivalent  Fresn  1  reflection  coefficient 
for  reflection  from  a  tilted  planar  surface. 

Subroutine  WGTAMP  modifies  the  specular  reflection  parameters  by  incor¬ 
porating  the  effects  of  the  transmitting  antenna  and  uses  the  subroutine  GXMTR 
which  provides  the  model  for  the  transmitting  antenna.  SETVAL  transfers  the 
multipath  quantities  to  the  labelled  common  areas  MLFLTI  and  MULTIP.  Finally, 
the  subroutine  DIGLOD  transfers  the  maximum  amplitude  building  multipath  to 
the  labelled  common  area  DIAGP  after  it  has  been  weighted  with  the  receive  an¬ 
tenna  pattern  modelled  by  subroutine  GAIN. 

It  is  to  be  noted  that  for  the  tactical  scenarios  associated  with  the  EHT 
studies,  the  number  of  buildings  specified  in  BDlfHFM  was  normally  set  to  zero 
so  that  building  reflection  segment  was  not  invoked. 

Appendix  B  supplies  additional  information  on  the  subroutines  shown  in 
Fig.  3.7  and  discussed  briefly  above. 

3.4  Aircraft  Reflection  Computation 

Figure  3.8  provides  the  program  sequence  for  the  computation  of  aircraft 
reflections.  Comparing  with  Fig.  3.7  a  general  similarity  in  the  program  flow 
is  noted.  As  for  the  case  of  buildings,  four  different  bounce  paths  are  con¬ 
sidered  are  also  considered  for  aircraft.  For  the  aircraft  case,  however, 
separate  consideration  is  given  for  fuselage  components  and  vertical  tail  com¬ 
ponents.  Both  are  modelled  as  cylinders:  the  fuselage  as  a  horizontal  cylin¬ 
der  and  the  tail  as  a  vertical  cylinder  and  the  same  basic  programs  are  used. 
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Fig.  3.8.  Program  processing  for  aircraft  reflections. 
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Referring  co  Fig.  3.8,  Che  subroutine  CHKMLT  is  used  twice  to  establish 
that  the  number  of  multipath  components  does  not  exceed  a  specified  value. 
The  first  time  is  for  the  fuselage  components  and  the  second  time  is  for  the 
fin  components.  After  the  first  call,  the  subroutine  FREFC  is  called  to  make 
the  calculation  for  the  fuselage  components.  FREFC,  in  turn,  calls  a  number 
of  other  subroutines: 

ACTYPE  determines  the  geometrical  parameters  associated  with 
each  aircraft 

CYLREF  computes  the  multipath  parameters  of  a  vertical  cylin¬ 
drical  surface 

ORIENT  is  used  to  compute  the  equivalent  vertical  cylinder  so 
that  CYLREF  can  be  used  for  horizontal  cylinders. 

The  subroutine  CYLREF  calls,  in  turn,  a  number  of  other  subroutine': 

VSUB,  XDIV  these  subroutines  are  vector  manipulation  subroutines 

NORMLZ,  ROTATE 

ORIEND 

CSPR  as  noted  previously,  is  used  to  compute  the  Fresnel  in¬ 

tegral. 

POLZN  computes  the  equivalent  Fresnel  reflection  coefficients 
for  reflection  from  a  tilted  planar  surface. 

After  the  second  call  of  CHKMLT,  the  tail  multipath  contributions  are 
computed  using  TREFC.  Since  the  fins  are  vertical  the  subroutine  ORIENT  is 
not  called,  and  the  computations  controlled  by  TREFC  for  the  tail  are  similar 
to  those  for  the  fuselage. 

As  for  the  building  computations,  there  are  calls  to  WGTAMP,  SETVAL,  and 
DIGLOD.  The  subroutine  WGTAMP  modifies  the  specular  reflection  parameters  by 
Including  the  effects  of  the  transmitting  antenna  and  uses  the  subroutine 
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GXMTK  to  provide  Che  model  of  the  antenna  pattern.  SETVAL  transfers  the  mul¬ 
tipath  quantities  to  the  labelled  common  areas  MULTI  and  MULT1P.  Toe  subrou¬ 
tine  DIGLOD  writes  the  maximum  amplitude  building  multipath  to  the  labelled 
common  area  DIAGP  after  It  has  been  weighted  with  the  receive  antenna  pattern 
modelled  by  subroutine  GAIN. 

The  subroutine  MAXMLT  Is  used  to  determine  the  scaling  parameters  for 
multipath  diagnostics  at  a  particular  evaluation  point.  A  call  to  PT1ME  pro¬ 
vides  the  time  marker  needed  to  compute  the  elapsed  time  employed  by  the  total 
time  needed  to  confute  the  multipath. 

finally,  the  subroutine  UIGOUT  writes  the  multipath  data  to  the  disk. 

As  for  the  case  of  building  reflections,  tactical  scenarios  associated 
with  KHT  studies,  did  not  normally  invoke  reflections  from  aircraft  so  the 
number  of  aircraft  in  BDUHFM  was  set  to  zero. 

Appendix  B  supplies  additional  data  on  the  subroutines  shown  in  Fig.  3.8. 

5.5  Parameters  Required  for  Multipath  Computations 

The  following  parameters  specify  the  scenario  which  is  used  in  the 
modelling  of  the  signal  in  space  seen  by  the  EHT  seeker.  A  standard 
rectangular  coordinate  system  is  used.  All  lengths,  frequencies,  and  times 
are  given  in  feet,  Hz,  and  seconds,  respectively.  The  symbols  used  are 
consistent  with  those  given  in  ( 1 ) .  In  most  cases,  parameters  are  set  to 
their  default  values  so  as  to  decrease  the  complexity  of  the  multipath 
confutations. 

Transmitter  Parameters  (Azimuth,  DME,  Elevation  1,  and  Elevation  2) 

1.  X .Y.Z-coordinates  of  location  of  transmitter.  XMTRAZ 

2.  Wavelength  (nominally  0.2  feet  for  C-band).  WLAZ 

3.  Polarization  (vertical  or  horizontal).  PLZAZ  where 
TRUE  implies  vertical  polarization 

4.  Dimension  of  transmitter  antenna.  DEFAULT  ■  0  (except 
for  DIMDME  which  must  be  specified  as  nonzero) 
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5.  Number  of  intervals  when  transmitter  is  off,  which  is 

less  than  or  equal  to  10  NAZBLK.  Start  and  end  time  of 
blinking  interval  in  seconds.  TAZBLK.  (  ,  ) 

6.  Number  of  transmitters.  NXMTR 

7.  Relative  power  levels  of  emitters.  BXMPWR 

0.  Tracking  option  parameter  for  each  transmitter. 

KXMTRK  (  ) 

9.  Planar  azimuth  pointing  angle  for  each  transmitter  for 

tracking  option  2.  BAZTRK  (  ) 

10.  Planar  elevation  pointing  angle  for  each  transmitter 

for  tracking  option  2.  BELTRK  (  ) 

Specular  Ground  Reflection  in  Labelled  Common  SURFAC 

1.  Number  of  rectangular  surface  elements,  which  must  be 
less  than  or  equal  to  10.  NR 

2.  X,Y,Z-coordinate  of  two  corners,  plus  X,Z-coordlnates 

of  a  third  corner,  in  increasing  order  of  magnitude  for 
the  X-coordinate,  for  each  rectangular  surface  ele¬ 
ment.  SURFS  (  ,  ) 

3.  The  real  and  imaginary  relative  dielectric  constants, 

denoted  by  Eg,  -tj,  respectively,  and  the  root- 
mean-square  roughness  height,  o^,  for  each  rectan¬ 
gular  surface  element.  ERS  (  ,  )  SH2S  (  ) 

4.  An  indicator  parameter,  for  each  rectangular  surface 
element,  which  is  equal  to  0  or  1  depending  on  whether 
its  contribution  to  the  specular  ground  reflection  is 
to  be  considered  a  primary,  or  secondary,  component, 
respectively  [1].  The  total  number  of  such  components 
is  specified  to  be  less  than  or  equal  to  ten.  NRSPEC 
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5.  Number  of  triangular  surface  elements,  which  is  leas 

than  or  equal  to  10,  NT 

6.  X ,Y, Z-coordinates  of  the  three  corners  of  each  trian¬ 

gular  surface  element,  in  Increasing  order  of  magnitude 
for  the  X-coordinate.  SURFS  (  ,  ) 

7.  tR,  “ei,  °h»  for  each  triangular  surface  ele¬ 
ment.  ERS,  SH2S 

8.  An  indicator  parameter,  for  each  triangular  surface  el¬ 
ement,  which  is  equal  to  0  or  1  depending  on  whether 
its  contribution  to  the  specular  ground  reflection  is 
to  be  considered  a  primary,  or  secondary,  component, 
respectively  [1],  NTSPEC 

9.  Default  values  of  eR,  -ej,  oh,  which  are  used  in 

those  regions  not  specified  by  the  previously  defined 
rectangular  and  triangular  areas.  ERO,  SH20 

10.  Number  of  Fresnel  zones  used  in  integration  region, 

nominally  equal  to  2.8.  NFZ 

11.  Number  of  grid  points  used  along  the  major  and  minor 
axes  of  the  integration  regions,  which  are  equal  to  25 
and  11,  respectively.  NA  and  NB 

12.  A  parameter  which  is  0  if  the  full  numerical  calcula¬ 
tion  of  the  Fresnel  integral  is  to  be  calculated,  and 
is  1  if  the  simplified  Fresnel  reflection-roughness  co¬ 
efficient  is  to  be  computed  using  FRERFS.  ISPGRD 
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Scattering  from  Buildings  in  Labelled  Common  BUILD  BDGRND 

1.  Number  of  buildings  which  should  be  less  than  or  equal 
to  10.  NBLD 

2.  Height  of  bottom  edge  of  front  face  of  building  above 

ground,  H5otCom,  for  each  building.  HBOT  (  ) 

3.  Height  of  building,  relative  to  bottom  edge,  Kg,  for 

each  building.  HB(  ) 

4.  X,Y-coordinates  of  left-hand,  and  right-hand,  edge  of 

front  face  of  building,  for  each  building.  BLD  (  ,  ) 

5.  eg,  -cj.  Oh  for  each  building.  DICSTBC),  SH2B(  ) 

6.  eR,  -ej,  oh  for  ground  reflection. 

7.  Tilt  angle  of  building.  This  angle  is  positive  if 

building  is  above,  and  tilts  away  from,  the  centerline, 
or  if  building  is  below,  and  tilts  toward,  the  center- 
line.  TILT  (  ) 

8.  Differential  height  factor  of  ground  (the  height  of 

ground).  This  factor  is  positive  if  the  ground  is 
above  zero  height  level,  and  negative  otherwise. 

GRNDBD  (  ) 


Scattering  From  Aircraft  ACS 

L.  Number  of  aircraft,  which  should  be  less  than  or  equal 
to  10.  NAC 


i 
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.  Aircraft  type,  for  each  aircraft,  e.g.,  1(»  747), 

2<-707-320B),  3(-  727),  4<-  DC10),  5<«  C-124),  6(- 

Convalr  880),  7(-  Hastings),  A  subroutine  ACTYPE  Is 
called  using  the  appropriate  aircraft  type.  This  sub¬ 
routine  Is  contained  in  program  JCHOJ.  It  loads  the 
following  aircraft  related  quantities  Into  storage. 
NACTYP  (  ) 

2.1  Area  of  both  wings.  WA 

2.2  Radius  of  fuselage.  FR 

2.3  Length  of  fuselage.  FL 

2.4  Radius  of  curvature  of  tail  fin.  TR 

2.5  Width  of  tail  fin.  TL 

2.6  Height  of  tall  fin.  TH 

2.7  Height  of  center  of  fuselage  above  the  ground.  HF 

2.8  Length  of  one  wing.  WL 

2.9  Width  of  wing.  WW 

2.10  Thickness  of  wing.  WT 

Quantities  2.8  through  2.9  are  found  In  labelled  common  WINGS. 

3.  X  ,Y-coordinates  of  cockpit  and  tall  fin  edge  of  fuse¬ 
lage  centerline,  for  each  aircraft.  AC  (  ,  )AC 

4.  Altitude  of  aircraft  for  each  aircraft,  defined  as  the 
height  of  fuselage  centerline  above  the  ground. 

5.  Differential  height  factor  of  the  ground. 

Shadowing  Due  to  Runway  Hump  in  LABELLED  COMMON  HMPR 
1.  X.Y.Z-coordinates  of  front  of  hump.  HUMPF 
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X,Y,Z-coordinates  of  center  of  hump.  HUMP 


3.  X,Y ,Z-coordlnates  of  back  of  hump.  HUMPB 

4.  Parameter  which  is  zero  if  runway  hump  shadowing  is  not 
to  be  done  and  one  otherwise.  IHUMP 

The  following  information  is  also  specified: 

1.  Length  of  runway.  RUNLEN 


2. 

Width  of 

runway. 

RUNWID 

Default 

Specif icati 

on  of  Multipath  Quantities 

1. 

DIMDME 

Nonzero 

2. 

ISPGRD 

1 

compute  specular  reflection  coefficients  using 

fast  confutation. 

3. 

NR 

0 

no  rectangles 

4. 

NT 

0 

no  triangles 

5. 

NB 

0 

no  buildings 

6. 

NAC 

0 

no  aircraft 

7. 

IHUMP 

0 

no  humps  and  shadowing  effects 
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4.0  SEEKER  MODULE 

This  section  will  discuss  In  a  general  way  the  program  flow  of  the  seeker 
module  as  was  done  In  Sections  2.0  an  3.0  for  the  Initialization  and  signal  In 
space  modules,  respectively.  Three  main  elements  make  up  the  seeker  module: 
(1)  the  antenna  model,  (2)  the  signal  processing  algorithm  and  (3)  the  tracker 
algorithm.  In  addition  to  the  description  of  the  program  flow,  subsections 
4.2  and  4.3  will  provide  some  additional  Insight  Into  the  signal  processing 
and  tracker,  respectively.  More  detail  on  the  various  subroutines  will  be 
given  In  Appendix  C,  along  with  definitions  of  the  variables  used  In  the  pro¬ 
gram  listings  which  accompany  this  report. 

4.1  Program  Sequence  for  Seeker  Module 

The  program  sequence  for  the  seeker  module  Is  shown  In  Fig.  4.1.  As  pre¬ 
viously  noted,  the  call  to  this  module  is  from  the  master  simulation  driver 
SJCMLST.  The  seeker  module  models  the  received  antenna  characteristics,  sig¬ 
nal  processing  and  angle  filtering.  The  input  to  this  module  include  the  am¬ 
plitude,  phase,  time  delay,  planar  incident  angles,  and  Doppler  frequency 
shift  of  each  of  the  multipath  components.  The  net  output  from  the  seeker 
module  is  two  planar  angles  in  azimuth  and  elevation  which  describe  the  orien¬ 
tation  of  the  line-of-sight  (LOS)  of  the  vehicle  to  the  estimated  location  of 
the  target.  These  are  used  in  the  guidance  computations. 

The  simulation  driver  calls  HOJ  which  in  turn  calls  the  subroutines  which 
model  the  signal  processing.  Subroutine  HOJ  in  turn  calls  RDPHAS  and  DLS  if 
specified  in  the  EXEC.  Subroutine  RDPHAS  adds  a  random  phase  between  -  and  2w 
to  the  phase  of  each  of  the  nultipath  components.  Subroutine  DLS  calls  C1HOJ 
which  models  the  enhanced  interferometer  processor.  This  DF  estimator  gener¬ 
ates  unbiased  estimates  on  the  strongest  emitter  in  the  f ield-of-view  by  per¬ 
forming  an  average  of  the  phase  difference  between  two  pairs  of  antenna  ele¬ 
ments  which  are  used  to  find  elevation  and  azimuth  conical  angles.  Note  that 
the  subroutine  GAIN  called  by  C1H0J  models  the  F.HT  antenna  characteristics 
given  the  angle  of  arrival  of  the  multipath  signals.  The  details  of  the  C1H0J 
subroutine  are  discussed  in  Subsection  4.2  below.  The  outputs  of  the  signal 
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Fig.  4.1.  Seeker  program  sequence. 
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processing  block  are  planar  angles  which  are  corrupted  by  errors  which  are 
characterized  in  ERRMOD.  This  subroutine  models  DF  angle  biases,  antenna 
measurement  errors  and  blind  range  effects.  If  specified,  the  corrupted  pla¬ 
nar  angle  estimates  are  filtered  by  the  tracker  algorithm  in  subroutine 
TRACK.  Further  description  of  the  tracker  algorithm  is  given  in  subsection 
4.3  below.  Finally,  subroutine  PUTHOJ  is  called  to  write  out  the  vehicle  po¬ 
sition,  velocity,  boresight  coordinates,  planar  angle  estimates  from  the 
tracker,  error  angles,  and  planar  difference  angles  to  the  disk. 

4.2  Enhanced  Interferometer  Signal  Processing 

The  signal  processor  is  based  on  the  phase  interferometer  concept.  The 
module  computes  the  real  and  imaginary  parts  of  each  signal  present  in  the 
electromagnetic  environment,  at  each  of  the  antenna  elements.  The  number  of 
signals  that  are  modelled  is  dependent  on  the  number  of  emitters  present  and 
the  number  of  multipaths  associated  with  each  emitter  using  amplitude  and 
phase  data.  The  phase  is  calculated  using  quantities  characterizing  multipath 
phases,  time  delays,  fractional  doppler  frequencies,  and  the  changes  in  phase 
due  to  differential  path  length.  Changes  in  phase  due  to  the  characteristics 
of  the  receive  antenna  can  also  be  incorporated.  The  real  and  imaginary  parts 
of  each  signal  sensed  by  the  antenna  elements  are  calculated  using  the  phase 
and  amplitude  of  each  multipath  conponent,  and  the  gain  characteristics  of  the 
receive  antenna.  Antenna  characteristics  are  modelled  by  calling  subroutine 
GAIN.  This  call  is  only  implemented  if  IGAIN  is  1  in  BDUHFM.  The  gain  is,  of 
irse,  dependent  on  the  angle  of  arrival  of  each  multipath.  The  processing 
scheme  employs  averaging  the  phase  difference  over  the  30  samples  taken  in  0.1 
sec.  Therefore,  signal  characteristics  are  determined  for  each  sample  in  the 
100  ms  period  by  Incorporating  the  appropriate  doppler  effects,  and  leaving 
all  other  multipath  quantities  unchanged.  It  should  be  noted  that  by  modify¬ 
ing  the  overall  simulation  '.nterval  (quantity  DRATE  in  block  data  file 
BDUHFM),  and  the  number  of  phases  generated  per  call  (quantity  FP  in  C1H0JG2), 
the  realism  of  the  model  is  increased.  This  occurs  because  the  effect  of 
changes  in  vehicle  position  and  orientation  on  multipath  attributes  over  the 
interval  that  the  50  signal  samples  are  collected  Is  better  simulated. 
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The  effects  of  measurement  noise  on  the  received  signal  In  space  can  also 

v 

be  modelled.  The  noise  has  a  Gaussian  distribution.  The  magnitude  of  the 
noise  Is  dependent  on  the  distance  of  Che  vehicle  from  the  emitter  specified 
as  No.  1. 

The  phase  difference  (for  each  of  the  50  samples  considered)  between  the 
azinuth  (1  and  3)  and  elevation  (2  and  4)  elements  Is  next  calculated.  The 
phase  difference  for  the  azimuth  channel  Is  resolved  about  the  quantity  AMI  5 
(specified  in  block  data  file  BDUHFM)  and  u .  The  phase  difference  between 
elements  2  and  3  Is  next  determined. 

The  phase  differences  between  elements  (2,  3)  and  (2,  4)  are  used  in  com¬ 
bination  to  resolve  ambiguities  in  the  estimate  of  the  phase  difference  be¬ 
tween  elements  (1,  3).  Phase  differences  are  averaged  over  for  azimuth  and 
elevation  channels.  Phase  differences  are  then  resolved  about  the  averages 
and  n  to  remove  ambiguities  and  then  averaged  again.  Transformation  to  coni¬ 
cal  angles  in  pitch  and  yaw  is  performed  using  the  known  wavelength  of  the 
emitters  and  the  antenna  element  separation.  (This  assumes  perfect  calibra¬ 
tion  of  the  antenna.)  The  conical  angles  are  then  converted  to  planar  angles 
which  are  used  In  the  models  for  tracker  and  guidance.  Angle  estimates  are 
modified  In  a  call  to  an  error  modelling  subroutine  prior  to  an  optional  call 
to  the  tracker  (pre-filter).  A  block  diagram  of  the  processing  Is  given  in 
Fig.  4.2  for  further  clarity.  j 

4.3  Tracker  Algorithm 

The  tracker  module  serves  to  filter  the  planar  angle  estimates  generated 
by  the  signal  processor  prior  to  their  use  in  guidance  software.  The  software 
performs  a  number  of  different  functions.  It  compares  the  D.F.  estimates  an¬ 
gles  vs.  truth  angles  for  each  of  the  emitters  that  exist  in  the  field.  Pla¬ 
nar  angles  which  are  body  referenced  (since  the  EHT  seeker  Is  body  fixed)  are 
transformed  using  a  detailed  Euler  transformation  matrix  (or  simpler  methods) 
to  space  referenced  angles.  A  number  (N)  of  the  space  referenced  elevation 
and  azimuth  angles  are  accumulated.  The  collected  data  Is  used  to  calculate 
'best *  tracker  estimates  of  the  target  position  using  a  1  coefficient  least 
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ENTER 


Fig.  4.2  Signal  processing  sequence. 
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squares  fit.  This  fit  is  made  either  through  a  brute  force  calculation  or 
through  an  iterative  technique.  Subsequent  estimates  (in  both  pitch  and  yaw) 
are  declared  good  if  the  space  referenced  azimuth  estimate  passes  a  variance 
test  which  uses  the  magnitude  of  the  variance  and  the  'best'  tracker  estimate 
in  azimuth.  If  the  estimate  passes  this  test,  it  is  used  (together  with  the 
N-l  previous  'good'  estimates)  to  rederive  best  tracker  estimates  and  the  var¬ 
iance  magnitude.  The  original  -  itch  and  yaw  body  referenced  estimates  are 
passed  to  guidance  without  modification.  If  the  azimith  angle  fails  the  vari¬ 
ance  test,  it  (together  with  the  corresponding  elevation  estimate)  is  replaced 
by  angles  computed  by  using  the  best  tracker  estimates  and  transforming  back 
into  the  body  frame.  The  tracker  then  enters  coast  mode.  Coast  logic  depends 
on  the  state  of  the  vehicle  (altitude  hold  or  pitchover  initiated). 

A.  Altitude  Hold 

For  a  finite  period,  AT,  all  estimates  will  be  rejected  and  replaced  by 
tracker  generated  angles.  Fallowing  this  interval,  if  N  successive  estimates 
pass  a  footprint  test,  and  the  variance  of  this  data  is  small  relative  to  a 
prespecified  constant  a,  the  vehicle  will  reenter  track  mode  and  estimates 
will  be  passed  unmodified  to  the  guidance  computer.  (The  footprint  test  de¬ 
termines  whether  the  vehicle  can  home  on  the  target  corresponding  to  the  new 
estimates. ) 

B.  Pitchover  Initiated 

If  the  vehicle  is  in  the  pitchover  mode,  the  rejection  of  one  set  (Az, 
El)  of  estimates,  does  not  automatically  entail  the  rejection  of  subsequent 
estimates.  'Best'  tracker  estimates  and  the  variance  are  frozen  until  N  suc¬ 
cessive  angles  pass  the  variance  criterion. 


5.0  FLIGHT  DYNAMICS  MODULE 

The  program  sequence  for  the  flight  dynamics  module  is  shown  in  Fig. 
5.1.  As  noted  in  Section  1.0  the  call  to  the  flight  dynamics  module,  FLTDYN, 
is  made  from  the  master  simulation  driver  $JCMLST.  FLTDYN  controls  the  se¬ 
quencing  of  the  other  subroutines  in  the  module: 

CRSGID  Guidance  subroutine 

AUTPLT  Autopilot  subroutine 

AERDYN  Aerodynamics  subroutine 

EQAMD  Equations  of  motion  subroutine 

KINEM  Kinematics  subroutine 

As  may  be  seen  in  the  diagram  some  of  these  subroutines  call  still  other  sub¬ 
routines.  This  section  will  discuss  in  a  general  way  the  purpose  of  these 
subroutines,  while  Appendix  D  will  provide  additional  details  of  the  input  and 
output,  as  well  as  the  definition  of  the  variables  as  given  in  the  program 
listings  which  accompany  this  report. 

One  of  the  first  tasks  of  the  FLTDYN  subroutine  is  to  initialize  the 
parameters  in  the  subroutines  which  it  calls.  In  particular,  it  calls 
separately  the  subroutine  WIND,  which  models  the  wind  shear  and  turbulence  to 

which  the  vehicle  is  exposed.  FLTDYN  also  controls  the  cycling  interval  of 

the  various  subroutines.  For  example,  the  cycling  rate  for  the  guidance 

computer  (CRSGID)  is  slower  than  for  the  other  subroutines  called.  In 

addition,  it  is  to  be  noted  that  the  cycling  rate  for  the  FLTDYN  subroutine 
will,  in  general,  be  different  than  for  the  subroutines  of  the  multipath  part 
of  the  simulation  discussed  in  Section  3.0. 

After  initialization,  the  first  subroutine  to  be  called  is  CRSGID,  the 
guidance  conq>uter  module.  The  key  inputs  to  the  guidance  computer  are  the 
planar  angle  measurements  which  specify  the  orientation  of  the  vehicle  with 
respect  to  the  estimated  line-of-sight  (LOS)  to  the  target  and  the  pitch  and 
yaw  body  angular  rates.  These  rates  are  filtered  and  combined  to  determine 
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Fig.  5.1.  Flight  dynamics  program  processing. 
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pitch  and  yaw  commands  to  the  autopilot  using  both  pursuit  and  proportional 
navigation  laws.  Filtered  planar  angles  are  also  used  as  a  basis  for  specify¬ 
ing  mode  switch  conditions  to  the  autopilot  to  transition  from  altitude  hold 
to  dive  on  the  simulated  target.  Several  seconds  after  pitchover  initiation, 
a  guidance  mode  switch  occurs  and  the  vehicle  transitions  from  pursuit  to  pro¬ 
portional  navigation.  The  guidance  computer  also  calls  the  subroutine  MICRO 
which  simulates  the  effects  of  the  actual  word  length  used  in  the  micro¬ 
computer  on  which  the  actual  guidance  computer  is  implemented. 

The  next  subroutine  called  by  FLTDYN  is  the  autopilot  module.  The  key 
inputs  to  the  autopilot  module  are  the  pitch  and  yaw  navigation  commands,  the 
mode  switch  commands  and  the  pitch  and  yaw  body  rates.  The  commands  are 
provided  by  the  guidance  module  while  the  pitch  and  yaw  body  angular  rates  are 
confuted  in  the  equation-of-motion  module  (EQAMO).  The  outputs  are  the  con¬ 
trol  surface  deflections  which  specify  the  orientation  of  the  elevator,  rudder 
and  ailerons.  It  is  to  be  noted  that  the  pitch  and  yaw  body  rates  provided  to 
both  guidance  computer  and  autopilot  are  corrupted  with  errors  which  model  the 
actual  gyro  performance. 

After  return  from  the  autopilot  module,  FLTDYN  calls  the  aerodynamics 
module,  AERDYN,  which  computes  the  aerodynamic  coefficients  and  thrust  on  the 
vehicle.  The  aerodynamic  coefficients  are  nonlinear  functions  of  many  varia¬ 
bles  including  angles -of -attack  and  sideslip,  control  surface  deflection,  body 
rates,  and  thrust.  In  addition,  the  linear  and  angular  velocities  associated 
with  the  wind  field,  which  are  computed  in  the  equations-of -motion  module,  are 
also  input.  It  is  further  noted  that  the  module  also  linearly  increases  the 
RPM  of  the  engine  after  pitch  over  to  reflect  the  observed  engine  speed  up 
that  occurs  in  a  dive.  Since  the  characteristics  of  the  thrust  are  given  in  a 
two-dimen8tional  table  (RPM,  velocity)  a  two-dimensional  table  look-up  rou¬ 
tine,  ATRPL2 ,  is  called  by  AERDYN. 

After  returning  from  the  aerodynamics  module  FLTDYN  calls  the  equation- 
of-motion  module,  EQAMO.  This  subroutine  uses  the  aerodymamic  coefficients 
produced  by  AERDYN  and  the  inertial  parameters  of  the  vehicle  to  calculate  the 
linear  and  angular  accelerations.  The  angular  accelerations  are  integrated  to 
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provide  angular  rates  (using  INTGRT)  which  are  then  used  to  compute  the  Euler 
or  space  angular  rates.  The  Euler  rates  are  then  integrated  (again  using 
INTGRT)  to  give  the  Euler  angles  of  the  body  defining  its  orientation  in 
space.  Also  in  this  module  the  linear  and  angular  components  of  the  wind 
field  in  the  vehicle  coordinate  system  are  computed  from  the  subroutine  WIND. 
In  accomplishing  this  task,  WIND  calls  RCWIND  which  transforms  the  inertlally 
referenced  wind  field  into  body  reference  components.  In  turn  RCWIND  calls 
RCVANG  which  uses  inertlally  referenced  coordinates  to  determine  the  azinuth 
and  elevation  angles  with  respect  to  the  vehicle. 

The  BLIMIT  call  used  by  the  EQAMO  subroutine  is  used  to  limit  the  eleva¬ 
tion  Euler  angle  (0)  from  coming  too  close  to  90  degrees  (n/2)  since  in  the 
computation  of  Euler  angle  a  division  by  cos6  is  made. 

Finally,  a  subroutine  to  compute  gyro  errors  GYERR  is  called.  This  sub¬ 
routine  provides  the  measurement  of  pitch  and  yaw  body  rates  corrupted  by 
bias,  scale  factor,  and  minimum  detectable  signal  level  which  is  used  by  both 
the  guidance  computer  and  autopilot. 

The  last  call  made  by  FLTDYN  is  to  the  subroutine  KINEM  which  takes  the 
linear  velocities  in  the  vehicle  frame  and  transforms  them  to  linear  veloci¬ 
ties  in  the  inertial  frame  using  a  transformation  matrix  computed  from  the 
Euler  angles  (a  call  to  INVERT).  These  vehicle  space  velocities  are  inte¬ 
grated  (a  call  to  INTGRT)  to  give  vehicle  space  positions. 
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6.0  THE  USE  OF  THE  EMITTER  HOMING  SIMULATION 

The  emitter  homing  simulation  is  run  in  time-share  mode  by  typing  in  the 

EXEC, 

TESTSMM  EXEC  A 

This  EXEC,  together  with  the  EXEC  TESTSIM2,  loads  the  key  programs  for  the  ex¬ 
ecution  of  the  simulation.  The  EXEC  TESTSIM2  also  serves  to  specify  the  ini¬ 
tial  conditions  for  the  simulation  case  considered. 

Most  of  the  simulation  programs  do  not  model  the  tactical  scenario, 
seeker  processing,  and  vehicle  flight  dynamics.  Those  programs  that  charac¬ 
terize  these  parameters  and  critically  affect  the  output  of  the  simulation  are 
listed  below: 

(1)  BDUHFM 

(2)  GDDRONE 

(3)  GDAUT 

(4)  GDGUID 

(5)  CSRPVTRK. 

(6)  C1H0JG2 

(7)  CCGAIND 

Those  programs  required  to  run  the  EHT  simulation  are  shown  in  Table  6.1.  The 
outputs  from  the  simulation  program  for  a  typical  run  (no  wind  effects,  no  an¬ 
tenna  measurement  errors,  no  gyro  errors)  are  included  in  the  next  section. 
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TABLE  6.1 

terminal  activity  with  definitions 


.  axac  a 

E*C  TESTS  IM2  0DUHFM  GOORONE  GOAUT  GOGUIO  CATI  OSRPVTRK  C  IHOJG2  CCGAINO 

R;  T<.01/  .*1  12:22:57 

. t  axac  a 

4CONTROL  AU 

E*C  S 

EXEC  8L  IPOF 

•e*c  TOISK  D 

EXEC  F1QTTP  TAPI 
EJCC  FI  40  A 
EXEC  FI®  A 
E*C  F !  60  A 
EXEC  CSGAIN 

LOAD  IJCH.ST  41  42  45  44  45  46  V  48  CRSPLTI  CR5PLT2  HOJRDT  CLSLOOP  M.TPLOT  JCHOJ  ©SCOPES  PNT30  (NOlAP 
4BEGSTAOC 

(A)  1090  ENTER  RUN  ID  NIMBER 

(B>  «QM-106  RUN  553AHZ  8KSP  3EMITTERS  REAL  ESTIMATE  DATA 

(C)  YES  DO  YOU  WANT  SPECULAR  GROUND  REFLECTION  ?<YES:NO> 

(D)  TOG  TYPE  IN  TDG  TO  00  I  NTERFEROETER  PROCESSING 

<E)  ENO  TYPE  IN  ENO  TO  TERMINATE  RCVP  TYPE  I  tf*UT 
<F>  13300  -330  2330  type  in  initial  position,  ft 

<G>  150  -2  »2  type  in  SPEED, FT/SEC,MI$SAZ, MISSEL,  OEWEES 

IH>  TEXC 
IE  NO 
START 

FI  10  CLEAR 
FI  14  CLEAR 
f I  i 5  Clear 

R;  W.01/0.03  13:23:04 


FILE  DEFINITIONS: 


SJCMLST 

8DUMFM 

GOORONE 

GOAUT 

GOGUID 

CATI 

CSRPVTRK 

C1H0JG2 

CCGAINO 

CR9»LT1 
CRSPL  T2 
HOJRDT 
CLSLOOP 
ML  TP  LOT 
JCHOJ 
JCSOOPES 
PNT30 


ENT  SWAJLATlON  0RIV7R. 
aoa  DATA  FILE, 

FLIGHT  OTNAMICS  DRIVER,  AERODYNAMICS,  EONS.  OF  ACTION. 

AUTOPILOT  MOOULE. 

guidance  mooule. 

SUBROUTINE  TO  CALL  ENHANCED  I  NTERFERCMETER  MOOULE. 

TRACKER. 

ENHANCED  I  NTERFEROETER  MOOULE. 

RECEIVE  ANTENNA  MDOULE. 

WO® AM  FOR  PLOTTING  OF  FLIGHT  OYNAMICS  ANO  SEEKER  PARAMETERS. 

WO® AM  FOR  SITING  ONTO  ANO  REN) INC  OFF  DISK, 

SUBROUTINE  FOR  DETERMINATION  OF  WHEN  TO  STOP  THE  SIMULATION, 

WO®AM  FOR  WRITING  OUT  l*»UT  INFLATION  ON  TACTICAL  SCENARIO  TO  TERMINAL. 
WO® AM  FOR  PLOTTING  OUT  MULTIPATH  QUANTITIES. 

WO® AM  FOR  ACOELLING  OF  TRANSMIT  ANTENNA  WTTERN. 

WO®AM$  FOR  SETTING  UP  PLOTTING  WO® AMS. 


ENTRY  DEFINITIONS: 

(A)  RUN  IDENTIFICATION  NIMBER 
(8)  TITLE  OF  RUN 

(C>  SPECIFICATION  OF  WHETHER  TO  USE  SPECULAR  GROW©  REFLECTION 
(0)  SPECIFICATION  OF  TYPE  OF  SIGNAL  WOCESSING  TO  BE  l*>LE»€NTED 

(E }  SPECIFICATION  ON  WHETHER  TO  L  IMIT  SIMULATION  RUN  TO  ONLY  ONE  S  IGNAL  PROCESSING  OPTION 
IF)  INITIAL  CHICLE  POSITION 

(G)  INITIAL  VEHICLE  SPEED,  AZIMUTH  ANO  ELEVATION  OFFSET  ANGLES  WITH  RESPECT  TO  LINE 
OF  SIGHT  TO  EMITTER  1. 

(H)  TEXC  INOICATES  TYPE  OF  TERMINAL  BEING  USED 
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7.0  FINAL  OUTPUT 


When  Che  vehicle  has  either  impacted  or  the  simulation  time  limit  has 
been  reached,  final  miss  distances  are  written  to  the  terminal  and  diagnostic 
plots  are  generated.  A  sample  set  of  annotated  plots  is  given  in  this  sec¬ 
tion.  Since  graphics  programs  are  not  interchangeable  between  the  Lincoln  and 
other  systems,  the  plotting  programs  listed  below  are  not  discussed  in  detail: 


Subroutine  MPOUT  (A,B) 
Program  MPOUT 


This  subroutine  outputs  the  minimum  miss  distance  for  a  flight  path  for 
the  closed  loop  simulation. 


Subroutine  AKPTMP 
Program  CLSLOOP 


This  subroutine  plots  the  HOJ  tactical  scenario  and  gives  the  amplitude 
rankings  of  the  multipath  components. 


Subroutine  FLTPUT  (A,B) 
Program  CRSPLTI 


This  subroutine  writes  out  flight  diagnostic  parameters  such  as  guidance 
commands,  accelerations,  body  rates,  and  Euler  angles. 


Subroutine  PLDT3D 
Program  CRSPLTI 


This  subroutine  plots  missile  trajectories  in  three  dimensions  and  in  two 


dimensions 


This  subroutine  plots  out  multipath  diagnostics. 


Subroutine  DATPLT 
Program  CRSPLT1 


This  subroutine  plots  out  planar  angles,  planar  projections  of  the  vehi 
cle's  trajectory  and  velocity  components. 
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1  CR  TO  GO  ON. . 

PROGRAM  TO  DO  CLOSED  LOOP  SIMULATION  FOR  EHT  TACTICAL  SCENARIO 
1000  XBQM-106  RUN  350MHZ  BKSP  3EMITTERS  REAL  ESTIMATE  DATA 


PARAMETERS 

FOR 

AZ  SYSTEM  ARE: 

PARAMETER 

VALUE 

UNITS 

COMMENTS  FOR  CLARIFICATION 

XMTRAZX 

0.0 

FT 

X 

1 

XMTRAZY 

0.0 

FT 

Y 

>  LOCATION  OF  EMITTER  1 

XMTRAZZ 

0.40000E+02 

FT 

Z 

1 

AZVELX 

0.0 

FT/S 

X  j 

i 

AZVELY 

0 .42000E+04 

FT/S 

Y 

•  VELOCITY  OF  EMITTER  1  (IRREVELANT) 

AZVELZ 

0.0 

FT/S 

Z  ' 

i 

WLAZ 

0.28200E+01 

FT 

WAVELENGTH 

PLZAZ 

T 

VERTICAL  POLARIZATION 

DIMAZ 

0.0 

FT 

DIMENSIONS  OF  TRANSMIT  ANTENNA 

ITYPAZ 

1 

EMITTER  1 

PARAMETERS 

FOR 

DME  SYSTEM  ARE : 

PARAMETER 

VALUE 

UNITS 

XMTRDMX 

0.69280E+04 

FT 

X  i 

I 

XMTRDMY 

0.40000E+04 

FT 

Y 

>  LOCATION  OF  EMITTER  2 

XMTRDMZ 

0 .50000E+02 

FT 

Z  ' 

1 

WLD 

0.28100E+01 

FT 

WAVELENGTH 

PLZD 

T 

VERTICAL  POLARIZATION 

DIMDE 

0.40000E+01 

FT 

DIMENSIONS  OF  TRANSMIT  ANTENNA  (MUST  BE  NONZERO) 

ITYPD 

2 

EMITTER  2 

PARAMETERS 

FOR 

ELI  SYSTEM  ARE: 

PARAMETER 

VALUE 

UNITS 

XMTRE1X 

0.69280E+04 

FT 

X  1 

XMTRE1Y 

-0.40000E+04 

FT 

Y 

|  LOCATION  OF  EMITTER  3 

XMTRE1Z 

0.15000E+02 

FT 

z  ! 

E1VELX 

0.0 

FT/S 

*  1 

1 

E1VELY 

0.0 

FT/S 

Y 

►  VELOCITY  OF  EMITTER  3  (IRRELEVANT) 

El  VELZ 

0.84000E+04 

FT/S 

z  ! 

1 

WLEl 

0.28300E+01 

FT 

WAVELENGTH 

PLZE1 

T 

VERTICAL  POLARIZATION 

DIME LI 

0. 1 1820E+02 

FT 

DIMENSIONS  OF  TRANSMIT  ANTENNA  (IRRELEVANT) 

ITYPE1 

3 

EMITTER  3 

PARAMETERS 

FOR 

EU  SYSTEM  ARE: 

PARAMETER 

VALUE 

UNITS 

XMTRE2X 

-0. 90000E+04 

FT 

X 

XMTRE2Y 

0 .25000E+05 

FT 

y  ; 

|  LOCATION  OF  EMITTER  4 

XMTRE2Z 

0 . 1 5000E+02 

FT 

Z 

E2VELX 

0.0 

FT/S 

X  1 

| 

E2VELY 

0.0 

FT/S 

Y  | 

>  VELOCITY  OF  EMITTER  4 

E2VELZ 

0.84000E+04 

FT/S 

Z  . 

) 

WLE2 

0.26250E+01 

FT 

WAVELENGTH 

PLZE2 

T 

VERTICAL  POLARIZATION 

DIMEL2 

0.78720E+01 

FT 

DIMENSIONS  OF  TRANSMIT  ANTENNA  (IRRELEVANT) 

ITYPE2 

4 

EMITTER  4 

Fig.  7.1.  Transmitter  parameters  as  listed  by  program  with  added  comments. 
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INPUTS  FOR  MODELING  OF  ERROR  SOURCES' 


Inputs  for  modeling  of  error  sources. 


C«  TO  M  Ml. 


r#s 


Fig.  7.5.  J’osition  and  velocity  coordin.it  cs.  final  miss  distances 
time  used  to  perform  multipath,  receiver,  vehicle  computations. 
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AMP  :  AMPLITUDE  IN  DB  OF  MULTIPATH  COMPONENT 

TIME  :  TIME  ASSOCIATED  WITH  MULTIPATH  COMPONENT 

RDOP  :  RECEIVER  FRACTIONAL  DIPPLER  FREQUENCY  THAT  OCCURS  AT  MAX.  WEIGHTED  AMPLITUDE 
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Fig.  7.8.  Elevator  angle/angle  of  attack,  pitch  body  rate 
Euler  angle  0. 
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Fig.  7.11.  Additional  elevation  and  azimuth  errors,  Elevation  of f-boresight  angle 
input  to  guidance  computer,  azimuth  of f-boresight  angle  input  to  guidance  computer 
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HOJ  MSfUC  TMJECTMIES 


Fig.  7.14.  Three-dimensional  view  of  vehicle  trajectory. 


Fig.  7.15.  Two-dimensional  projection  of  vehicle  trajectory. 
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Fig.  7.20.  Planar  azimuth  and  elevation  difference  angle  between 
Vehicle  boresight  vector  Relative  to 

Vehicle  velocity  vector  vector  to 

Seeker  estimate  vector  various  targets 
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Appendix  A 

SUBROUTINES  FOR  INITIALIZATION  MODULE 

This  appendix  provides  additional  details  relative  to  the  subroutines 
used  by  the  Initialization  module.  It  Is  most  useful  when  used  with  the 
listings  included  with  this  report. 

The  following  subroutines  are  included  along  with  the  page  number  on 
which  they  are  found: 


SUBROUTINE 

Page 

PTIME 

A- 2 

MISINT 

A-2 

READIT 

A-3 

RDRI 

A- 3 

YESNO 

A- 4 

RDRC 

A- 4 

PUTCHR 

A- 4 

WRTYN 

A- 4 

HOJRDT 

A- 5 

GETRCV 

A-6 

RDRCV 

A- 7 

RDRXYZ 

A- 7 

INTPOS 

A- 7 

INVERT 

A-8 

TITLPG 

A— 8 

BXMWRT 

A- 9 

INIT 

A- 10 

SURFSW 

A-10 

VSUB 

A-l  1 

NORML 

A-l  1 

CROX 

A-l  1 

NORMLZ 

A-l  1 

ORIENT 

A- 12 

DOT 

A-l  2 

ACTYPE 

A- 12 

SCLHOJ 

A-l  3 

FLIGHT 

A-13 

XDIV 

A-l  3 

MINPOS 

A- 14 

EMCASE 

A-l  4 

UPDTMS 

A- 14 

MOVE 

A-l  4 

INITR 

A-l  5 

ORIENT 

A-l  5 

A-l 


INITXR 

A- 13 

VSUB 

A-ll 

NVORMC 

A- 16 

VNORM 

A-16 

Subroutine  PTIME  (A,B,C,D,E,F,G) 
Program  PTIME 


This  routine  obtains  timing  information.  It  is  written  in  Assembler  for 
the  IBM  370.  AH  the  arguments  correspond  to  values  which  are  returned. 

A  =»  Integer,  in  milliseconds,  which  contains  the 
amount  of  CPU  time  used  by  the  task 
B  ■  Integer  giving  the  hour  of  the  day 
C  “  Integer  giving  the  minute  of  the  hour 
D  »  Integer  giving  the  second  of  the  minute 
E  »  Integer  giving  the  month 
F  »  Integer  giving  the  day  of  the  month 
G  *  Integer  giving  the  year 


Subroutine  MISINT 
Program  MISINT 


This  routine  initializes  parameters  LSTEST,  LTHRUS  as  FALSE  in  labelled 
common  TSTMIS  for  initializations  in  the  flight  dynamics  simulation  of  the 
MAVERICK  missile.  This  is  included  in  this  discussion  for  generality. 


A- 2 


Subroutine  READIT  (6) 
Program  $JCMLST 


This  subroutine  reads  the  number  which  indexes  the  run  and  the  title 
which  labels  the  run.  These  quantities  are  specified  in  the  EXEC.  The  EXEC 
input  is  also  used  to  determine  whether  specular  ground  reflections  will  be 
modelled  and  initial  conditions  for  the  vehicle  state  vector  (orientation, 
position,  velocity).  The  following  subroutines  which  are  discussed  below,  are 
used  to  read  information  from  the  EXEC: 

Subroutine  RDRI 
Subroutine  YESNO 
Subroutine  WRTYN 
Subroutine  HOJRDT. 

It  should  be  noted  that  if  ground  reflections  are  not  to  be  modelled  variable 
IFGREF  is  0;  if  they  are  to  be  modelled  IFGREF  is  1.  The  value  of  this  param¬ 
eter  determines  whether  subroutines  modelling  specular  ground  reflections  are 
to  be  called.  IFGREF  is  a  parameter  in  the  subroutine  READIT  and  is  set  by 
input  from  the  EXEC  which  is  read  by  subroutine  YESNO. 


Subroutine  RDRI  (I,J) 
Program  RDRI 


This  subroutine  reads  the  number  used  as  a  label  to  index  the  simulation 
run.  It  is  system  dependent. 

I  =*  5  which  corresponds  to  a  read  from  the  terminal 

J  »  Index  number  input  by  the  user 

If  J  is  input  incorrectly,  the  subroutine  returns  REENTER  until  the  input  is 
correct. 
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Subroutine  YESNO  (TORF) 
Program  YESNO 


This  subroutine  reads  the  user  response  to  the  question  of  modelling 
specular  ground  reflections.  TORF  is  a  logical  variable  which  is  set  to  TRUE 
is  the  answer  YES  is  given  to  the  question  and  FALSE  if  the  answer  NO  is 
given.  Either  YES  or  NO  must  be  input  in  response  to  this  question  otherwise 
the  subroutine  will  return  the  response  REENTER.  The  program  calls  RDRC  for 
the  reading  of  the  input  data  from  the  EXEC. 


Subroutine  RDRC  (I, A) 
Program  RDRC 


This  subroutine  reads  A  (either  'YES'  or  ’NO’)  from  the  EXEC.  If  A  is 
neither  ’YES'  or  'NO',  the  subroutine  PUTCHR  is  called  to  output  'REENTER'  to 
the  terminal.  I  is  set  to  5. 


Subroutine  PUTCHR  (fl) 
Program  PUTCHR 


This  subroutine  writes  the  work  'REENTER'  to  the  terminal.  B 
'REENTER'. 


Subroutine  WRTYN  (6,  TORF) 
Program  WRTYN 


This  subroutine  will  write  ’YES'  to  the  terminal  (unit  6)  if  specular 
ground  reflections  are  to  be  modelled,  and  'NO'  if  they  are  not  modelled  de¬ 
pending  on  the  value  of  the  logical  variable  TORF. 
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Subroutine  HOJRDT  (6) 
Program  HOJRDT 


This  subroutine  reads  which  receivers  are  to  be  considered  for  the  simu¬ 
lation  run.  The  capability  exists  to  model  vehicle  performance  in  the  speci- 
i  fied  tactical  scenario,  with  a  wide  variety  of  candidate  antenna/receiver/ 

signal  processing  systems  for  direction  finding.  It  also  accepts  as  input  the 
initial  position  of  the  vehicle,  and  its  speed  and  orientation  angles  with 
respect  to  Emitter  1.  Given  the  orientation  of  the  vehicle  velocity  vector 
relative  to  this  emitter,  the  coordinates  of  the  vehicle  boresight  with  re¬ 
spect  to  an  inertial  coordinate  system  is  determined. 

OUTPUTS 

LABELLED  COMMON 

H0J1  PIClC(I)  is  a  logical  array  which  speci¬ 

fies  which  DF  seekers  are  to  be  mod¬ 
elled.  Parameter  I  can  vary  from  1  to 
5.  PICK(2)  ■*  TRUE  if  and  only  if  RDG 
is  specified  in  the  EXEC  as  a  receiver 
option,  and  corresponds  to  the  use  of 
the  Enhanced  Interferometer. 

LABELLED  COMMON 

H0J2  POSITIONS)  -  Missile  position  (ft) 

VELCTY(3)  =»  Missile  velocity  (ft/sec) 

B0REST(3)  =*  Missile  boresight  (ft) 

These  coordinates  are  defined  relative 
to  an  inertial  system  and  are  constantly 
updated  throughout  the  simulation  run. 
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LABELLED  COMMON 


SPAHOJ  SPEED  «  Real  missile  speed  (ft/sec) 

MISSAZ  -  Real  initial  azimuth  planar 
angle  (degs)  relative  to  Emit¬ 
ter  1  (denoted  by  XMTRAZ) 
MISSEL  “  Real  initial  elevation  planar 
angle  (degs)  relative  to  Emit¬ 
ter  1  (denoted  by  XMTRAZ) 

LABELLED  COMMON 

HOJA  INP0S(3)  =»  Initial  missile  position 

INTVEL(3)  *  Initial  missile  velocity 

INTB0R(3)  ■  Initial  missile  boresight 


Subroutine  GETRCV  (PICK,  NUMRCV ,  PRINT,  OUT,  IN) 
Program  HOJRDT 


This  subroutine  determines  which  elements  of  the  logical  array  PICK  are 
TRUE  depending  on  the  receiver  options  which  are  specified.  TDG  must  be  spec¬ 
ified  in  the  EXEC  for  the  enhanced  interferometer  to  be  modelled.  When  this 
character  string  is  read  by  subroutine  KDRCV,  the  second  element  in  the  logi¬ 
cal  array  PICK  is  set  to  TRUE.  Other  character  strings,  when  input  in  the 
EXEC,  will  cause  other  elements  of  the  array  PICK  to  be  set  to  TRUE.  If  only 
TDG  is  specified  all  of  the  elements  of  PICK  are  set  to  FALSE  except  the  sec¬ 
ond  element  [PICK(2)J  which  is  set  to  TRUE. 

INPUTS 

In  argument  list  NUMRCV  =  5 

PRINT  =  b 
OUT  =  6 
IN  5 
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OUTPUTS 

In  argument  list 


PICK(I) 
noted  above. 


I 


1,5  determined  as 


Subroutine  RDRCV  (IN,  OUT,  RCV,  EOF) 
Program  HOJRDR 


This  subroutine  reads  the  characters  in  the  EXEC  file. 

IN  -  5 

OUT  -  6 

RCV  ■  Character  string  specifying  various  receiver  options 
EOF  =*  Logical  variable  which  is  TRUE  when  character  string  END  is 
read  in  the  EXEC  file. 


Subroutine  RDRXYZ  (X,Y,Z) 
Program  RDRXYZ 


This  subroutine  reads  in  three  real  numbers,  which  specify  the  vehicle 
position.  It  calls  ERRSET,  a  system  subroutine  for  checking  input  data  for 
errors. 


Subroutine  INTPOS  (Speed,  Az,  El,  Pos,  Vel,  Borest) 
Program  GDDRONE 


This  subroutine  uses  initial  specification  of  quantities  Speed,  Az,  and 
El  to  determine  the  vehicle  boresight  vector  coordinates  with  respect  to  iner¬ 
tial  space.  This  subroutine  calls  INVERT. 


This  subroutine  uses  information  regarding  the  line-of -sight  vector  from 
missile  to  Emitter  1  (POSX)  and  knowledge  of  the  velocity  vector  orientation 
relative  to  the  line-of-sight  (VELP)  to  determine-  velocity  coordinates  rela¬ 
tive  to  inertial  space  (VEL),  where  POSX,  VELP,  and  VEL  are  arrays. 


Subroutine  TL'CLPG 
Program  CLSLOOP 


This  subroutine  writes  out  some  of  the  parameters  associated  with  the 
emitters,  such  as  position,  wavelength,  and  polarization  on  the  first  page  of 
output.  On  the  second  page  of  output,  parameters  associated  with  the  multi- 
path  model  are  written  out.  Por  typical  cases  in  which  the  number  of  rectan¬ 
gular  and  triangular  surface  elements,  the  number  of  buildings,  and  tbe  number 
of  aircraft  is  zero,  and  the  number  of  emitters  is  3,  the  following  quantities 
are  written  out: 

XMTR _ X 

XMTR _ Y 

XMTR _ Z 

_ VELX 

_ VELY 

_ VELZ 

WL _ 

PL _ 

F  -  Horizonta 


Position  of  emitters 
in  ft. 

Velocity  of  emitters 
in  ft.  (not  incorporated 
into  the  simulation  model) 
Wavelength  of  emitters 
Polarization  T  •  Vertical 


DIM 


Dimensions  of  transmit  antenna  (Irrelevant 
except  that  DIMDME  must  not  be  0) 


ISPGRD 


0  if  full  integration  is  to  be  done 
1  if  simple  Fresnel  reflection  roughness  co¬ 
efficient  is  to  be  confuted 


ERO 

default  complex  dielectric  constant 

SH20 

default  rms  roughness  height 

NA 

25 

Integer  number  of  grid  points  along  major  and 

minor  axis  of  ground 

NB 

1 1 

reflection  integration  ellipse 

Parameters  for  diffuse  ground  scattering  are  next  written  out.  (It  should  be 
noted  that  diffuse  scattering  is  not  modelled  in  this  iteration  of  $JCMLST.) 
Parameters  for  editing  the  various  multipaths  are  also  output  to  the  termi¬ 
nal.  These  are  also  not  used  in  the  iteration  of  $JCMLST  being  presented. 
The  quantities  which  are  written  out  in  TITLPG  are  all  specified  by  the  user 
in  the  block  data  file  BDUHFM. 

It  should  be  noted  that  some  of  the  calls  in  this  program  assume  a  par¬ 
ticular  graphics  package.  Subroi  tines  PAGADV,  ORTSPC ,  OHMP,  OSHBD,  and 
OSHDAC,  which  are  called  by  TITLPG,  also  reside  in  program  CLSLOOP. 


Pub routine  BXMWRT 
Program  BXMWRT 


This  subroutine  writes  out  to  the  terminal  blinking  emitter  parameters 
initial  conditions  of  the  vehicle  and  parameters  associated  with  the  seeker. 


(1)  For  each  emitter,  the  relative  power,  relative  voltage,  and 
parameters  specifying  the  antenna  pattern  are  written  out. 

(2)  The  receiver  types  that  are  being  simulated  are  written  out 
to  the  terminal.  Since  TDG  is  specified  in  the  EXEC,  the 
receiver  that  is  modelled  is  the  interferometer.  Next  the 
variables 
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FAD2 


|  (meters)  Azimuth 


FED2 

-  —  (meters)  Elevation 

F1A  - 

2fly  (dimensionless) 

Azinuth 

FIE  = 

2fry  (dimensionless) 

Elevation 

are  output 
I  (JAIN,  and 

as  well  as  the  parameters  DYNX,  AMI  5,  ITRK, 
BWT.  These  parameters  are  used  in  the  modelling 

of  Che  enhanced  interferometer.  Lastly,  initial  vehicle 
parameters  such  as  speed,  misalignment  angles  in  azimuth  and 
elevation  between  the  L.O.S.  to  Emitter  I  and  the  velocity 
vector,  and  initiaL  position  and  velocity  components  are 
written. 

O)  The  number  of  blinking  intervals  and  the  off-time  intervals 
for  each  emitter  is  written  out. 


Subroutine  IMIT 
Program  INIT 


This  subroutine  initializes  the  common  areas  of  the  signal  in  space  rou¬ 
tines  based  upon  the  input  contained  in  the  block  data  file  BLHJHFM. 


Subroutine  SURFSW  (A,B) 
Program  SURFSW 


This  subroutine  modifies,  if  necessary,  the  SURFS  array  which  contains 
the  coordinates  of  the  corners  of  the  rectangles  and  triangles  that  are  used 
to  represent  the  ground  plane.  Other  subroutines  that  are  called  expect  that 
the  corners  will  be  in  a  specified  order.  This  subroutine  checks  the  SURFS 
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array,  and,  it'  necessary,  modifies  it  to  insure  that  the  proper  order  is  main¬ 
tained.  Integers  A  and  B  indicate  the  locations  in  the  SURFS  array.  When  the 
number  of  rectangles  and  triangles  is  set  to  zero,  as  is  typically  the  case, 
this  routine  is  not  called. 


Subroutine  VSUB  (A,B,C) 
Program  VSUB 


This  subroutine  subtracts  vector  B  from  vector  A  to  give  vector  C,  where 
A,  B,  and  C  are  arrays  of  order  3. 


Subroutine  NORML  (A,B,C) 
Program  NORML 


This  subroutine  computes  the  cross  product  of  vectors  A  and  B  and  normal¬ 
izes  to  produce  the  vector  C. 


Subroutine  CROX  (A,B,C) 
Program  CROX 


This  subroutine  computes  the  cross  product  of  vectors  and  A  and  B  and 
puts  the  result  in  C. 


Subroutine  NoRMLZ  (A) 
Program  NORMLZ 


This  subroutine  normalizes  vector  A. 


A-ll 


Subroutine  ORIENT  (A,B,C) 
Program  ORIENT 


This  subroutine  computes  the  components  of  a  vector  relative  to  a  verti¬ 
cal  cylinder,  given  its  components  with  respect  to  the  Ac^  horizontal  alr- 
cratt  Luseiage.  liiis  Clous  forma  Lion  is  employed  so  that  the  same  subroutines 
which  model  reflections  off  a  horizontal  aircraft  fuselage  can  be  used  to 
model  multipath  off  a  vertical  tall  fin.  This  subroutine  is  typically  not 
called  for  tactical  scenarios  and  the  number  of  aircraft  is  set  to  zero  in 
block,  data. 


Function  DOT  (A,B) 

Program  DOT 

_ 

This  function  computes  the  dot  product  of  two  vectors  with  coordinates  A 
and  B  . 


Subroutine  ACTYPE  (A ,B ,C ,D ,E ,F ,G ,H) 
Program  JCHOJ 


This  subroutine  specifies  the  following  parameters  for  certain  aircraft 
A  =  Aircraft  type 
B  =  Radius  of  fuselage 
C  =  Length  of  fuselage 
D  =  Radius  of  curvature  of  tail  fin 
E  =  Length  of  tail  fin 
F  =  Height  of  tail  fin 
G  =  Total  are  of  both  wings 

H  =  Height  of  center  of  fuselage  above  the  ground. 

When  the  number  of  aircraft  is  set  to  zero  this  routine  is  not  called. 
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Subroutine  SCLHOJ 
Program  SCLHOJ 


This  subroutine  provides  the  scaling  parameters  used  by  the  plotting  pro¬ 
grams. 


Subroutine  FLIGHT  (STOP,  NEWSEG) 
Program  HOJRDT 


This  subroutine  is  used  to  determine  when  a  particular  receiver  looi  is 
to  be  ended  and  final  miss  distance  calculated.  The  vehicle  trajectory  itera¬ 
tion  is  ended  when  one  of  several  criteria  are  met: 

(a)  The  index  which  specifies  the  flight  time  thru  the  Eq . 

FLIGHT  TIME  =  (Index)  (Drate) 

exceeds  a  certain  value,  IND,  which  is  specified  in  this 
subroutine 

(b)  The  altitude  of  the  vehicle  is  less  than  0.1  ft. 

(c)  The  altitude  of  the  vehicle  is  less  than  the  height  of  a 
particular  emitter  and  the  distance  from  the  emitter  is  less 
than  50  ft. 

If  any  of  these  criteria  are  met,  logical  variable  STOP  is  set  to  .TRUE.. 
Logical  variable  NEWSEG  is  always  set  to  .TRUE.. 


Entry  XDIV  (A,B) 
Program  NORMLZ 


This  subroutine  normalizes  the  receiver  velocity  vector  by  the  speed  of 
light. 
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Subroutine  MINPOS  ( IXMTR) 
Program  MINPOS 


This  subroutine  determines  the  minimum  miss  distance  relative  to  the 
emitter  IXMTR. 


Subroutine  EMCASE  ( IXMTR,  TEMP) 
Program  CRSPLT1 


This  subroutine  determines  the  coordinates  of  the  vector  TEMP  from  the 
vehicle  to  emitter  IXMTR. 


Subroutine  UPDTMS  (X,  IXMTR,  TEMP) 
Program  UPDTMS 


This  subroutine  determines  Che  miss  distance,  X,  from  the  vehicle  to  the 
emitter  IXMTR  using  a  point  of  closest  approach  analysis. 


Subroutine  MOVE  (A,B) 
Program  MOVE 


Subroutine  MOVE  is  used  to  load  emitter  parameters  into  labelled  common 

XMTR. 

A  -  arrays  BLKA ,  BLKD,  BLK2 
B  =  arrays  BLKT  in  common  labelled  XMTR 
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Subroutine  INITR 
Program  INITR 


This  subroutine  computes  the  receiver  position  and  velocity  coordinates 
in  the  reference  frame  of  equivalent  vertical  aircraft  fuselages. 


Subroutine  ORIENT  (A,B,C) 
Program  ORIENT 


This  subroutine  computes  the  components  of  a  vector  relative  to  a  verti¬ 
cal  cylinder  given  its  components  with  respect  to  a  horizontal  cylindrical 
fuselage. 


Subroutine  INITKR 
Program  INITXR 


This  subroutine  computes  the  relative  vector  from  the  transmitter  to  the 
receiver,  the  wave  number,  and  the  wave  number  squared. 


i 


Function  VNORMC  (A,B,C) 
Program  NVCRMC 


This  function  finds  the  length  of  a  vector  with  coordinates  A,B,C 
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Function  VNORM  (A) 
Program  VNORM 


This  function  finds  the  length  of  a  vector  with  coordinates 

clement  array,  A. 


A-16 


Appendix  B 

SUBROUTINES  FOR  SIGNAL-IN-SPACE  MODULE 


This  appendix  provides  additional  details  relative  to  the  subroutines 
used  by  the  signal-in-space  module  and  complements  Section  3.0.  It  is  most 
useful  when  used  with  the  program  listings  which  accompany  this  report.  For 
convenience  in  using  this  appendix  with  Section  3.0,  it  is  divided  into  four 
parts  B.l,  B .2 ,  B.3,  and  B.A  which  correspond  with  Sections  3.1  through  3. A. 

B.l  Direct  Wave  Subroutines 

The  following  are  the  subroutines  related  to  the  direct  wave  computation 
(discussed  in  Section  3.1)  along  with  the  page  number  on  which  they  are  found. 


Subroutine 

CHKMLT 

DIRWV 

RCVANG 

VNORMC 

DOTC 

ATAN22 

SHDHMJ 

CIRCSH 

GZ 

HMPROT 

CSPR 


Page  No. 
B-l 
B-2 
B-3 
B-3 
B-3 
B-A 

b-a 

B-b 

B-7 

B-7 

B-7 


Subroutine  CHKMLT  (A) 
Program  HOJRDT 


This  subroutine  checks  the  number  of  multipath  components  (A)  and  if  A  is 
greater  than  1A7  stops  execution  of  the  program. 
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Subroutine  DIRWV 
Program  DIRWV 


This  subroutine  determines  the  multipath  parameters  associated  with  the 
direct  wave. 

INPUTS 


LABELLED 

COMMON 

RCVR  (3)  - 

Receiver  position  relative  to  the 

RCVER 

RVEL  (3)  = 

transmitter  (3-vector). 

Receiver  velocity  (3-vector)  in  iner¬ 
tial  coordinates. 

LABELLED 

COMMON 

XMTR  (3)  * 

Transrai  .ter  position  (3-vector). 

X VEL  (3)  = 

Transmitter  velocity  (3-vector). 

OUTPUTS 

LABELLED 

COMMON 

DOPPD 

Doppler  of  direct  wave 

DIRWAV 

DOPPDR 

Doppler  of  receiver  (used  in  signal 

processor  module  to  model  the  phase 

change  in  the  multipath  components 

over  small  increments  of  time). 

AZDIP 

Azimuth  planar  angle  specifying  the 

location  of  the  receiver  with  respect 

to  the  transmitter  in  the  inertial 

reference  frame. 

ELDIR 

Elevation  planar  angle  specifying  lo¬ 
cation  of  receiver  with  respect  to  the 

transmitter  in  the  inertial  reference 

f  rame. 

AZIDIR 

Azimuth  planar  angle  relative  to  the 

vehicle  boresight. 

ELID1R 

Elevation  planar  angle  relative  to  the 

vehicle  boresight. 


ELID1R 


Subroutine  RCVANG  (A,B ,C ,D,E ,F ,G ,H) 
Program  GDDRONE 


X 

Y 

Z 


A 
B 
C 

D.E.F 
G  = 

H  = 


This  subroutine  uses  inertially  referenced  coordinates  to  determine  azi¬ 
muth  and  elevation  angles  relative  to  the  vehicle.  The  arguments  are 

vector  coordinates  to  a 
point  in  space  in  inertial 
reference  frame, 
dummy  arguments 

Planar  azimuth  angle  relative  to  vehicle 
boresight 

Planar  elevation  angle  relative  to  vehicle 
boresight. 


Function  VNORMC  (A.B.C) 
Program  VNORMC 


This  function  computes  the  length  of  a  three-component  vector  with  A,B,C 
being  the  X,Y,Z  coordinates. 


Function  DOTC  (A,B,C,D,EfF) 
Program  DOTC 


This  function  computes  the  dot  product  of  two  three-component  vectors 
with  A,B,C  being  the  coordinates  of  one  vector  and  D,E,F  the  components  of  the 
other. 
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Function  ATAN22  (A,B) 
Program  ATAN22 


I. 

t 


This  function  confutes  the  arc  tangent  of  A/B  and  returns  a  value  between 
-rr  and  n . 


Subroutine  SHDHMJ 
Program  SHDHMJ 


This  subroutine  models  the  shadowing  effects  due  to  circular  cylindrical 
humps,  and  determines  the  associated  multipath  parameters.  These  parameters 
are  like  those  found  in  DLRWV.  Note  that  this  subroutine  is  called  only  if 
IHLTMP  is  set  to  1  in  block  data  file  BDUHFM.  The  default  value  is  zero. 

KEY  INPUTS 
LABELLED  COMMON 

HUMPR  HUMPF  (3)  *  X,Y,Z  coordinates  of  the  front  of  cir¬ 

cular  cylindrical  runway  hump. 

HUMP  (3)  *  X,Y,Z  coordinates  of  the  middle  of 

circular  cylindrical  runway  hump. 

HUMPB  (3)  ■  X,Y,Z  coordinates  of  the  back  of  cir¬ 

cular  cylindrical  runway  hump. 

These  coordinates  assume  the  coordinate  system  shown  in  Fig.  3.4.  With  this 
coordinate  convention,  the  coordinates  HUMPF(2),  HUMP(2),  and  HUMPb(2)  are  all 
zero. 


LABELLED  COMMON 

XMTER  XMTR  (3)  =  X,Y,Z  coordinates  of  the  transmitter 

location. 


B-4 


LABELLED  COMMON 

RCVER  RCVR  (3)  “  X,Y,Z  coordinates  of  the  receiver  lo¬ 

cation. 

RELRCR  (3)  *  X,Y,Z  coordinates  of  the  vector  from 
receiver  to  the  transmitteL. 


LABELLED 

COMMON 

AIRPRT 

RNWID 

- 

Runway  width 

LABELLED 

COMMON 

0IRWAV 

AZDIR 

— 

' 

Planar  angles  (Azimuth  and  Eleva- 

i  tion)  of  the  direct  wave  in  the 

ELDIR 

= 

transmitter  coordinate  system. 

DOPPDR 

- 

Doppler  of  receiver 

AZIDR 

= 

Planar  angles  of  the  direct  wave 

■  as  seen  from  the  receiver  in  the 

ELIDR 

3 

receiver  coordinate  system. 

KEY  OUTPUTS 

LABELLED  COMMON  TSHAD(l)  <— >  Ed/Eq  where  ED/£0  is  ratio  of 

diffracted  wave  to  the  free  space  sig¬ 
nal. 

AZSHAD(l)  <— >  AZDIR 
ELSHAD(l)  <— >  ELDIR 
AZISHAD(1)<— >  AZ1DR 
ELI SHAD ( 1)< — >  ELIDR 
RD0P(1)  <— >  DOPPDR 

The  quantities  TSHAD,  AZSHAD,  ELSHAD,  AZISHAD,  ELISHAD  are  analogous  to  the 
multipath  parameters  determined  in  the  call  to  DIRWV.  The  quantity  TSHAD  is  1 
for  the  direct  wave  by  definition. 
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LABELLED  COMMON 


SHCIRC 


CENTER(2) 


Distance  from  the  center  of  the  hump 


to  the  vector  from  the  receiver  to  the 


transmitte  r. 


This  subroutine  carries  out  the  detailed  calculation  of  the  shadowing  ef¬ 
fect  due  to  a  circular  cylindrical  hump. 

A  =  Knife  edge  parameter  (a) 

B  =  a2 

C  =  Correction  term  to  be  subtracted  from  the  knife 
edge  result. 

(GOO  eO«2) 

u 

D  »  Specifies  the  statement  number  in  the  calling 
subroutine  that  this  subroutine  will  return  to 
if  certain  criterion  are  met. 

INPUTS 

LABELLED  COMMON 

RCVER  RCVR  (3)  =  Receiver  position  3-vector  in  the  in¬ 

ertial  reference  frame. 

RELRCR  (3)  *  3-vector  from  the  receiver  to  the 
transmitter  in  the  inertial  reference 
f  rarne . 
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XMTER 


XMTR  (3) 


Transmitter  position  3-vector  in  the 
inertial  reference  frame. 

Wave  number 


r 


k 


Function  GZ(A) 
Program  SHDHMJ 


This  is  the  complex  function  used  to  compute  the  shadowing  effect  due  to 
the  circular  cylindrical  hump. 


Subroutine  HMPROT  (A,B,C,D) 
Program  SHDHMJ 


This  subroutine  transforms  inertially  referenced  coordinates  of  a  vector 
into  a  coordinate  system  based  on  the  vector  from  the  transmitter  to  the  re¬ 
ceiver. 

A  -  cosine  of  angle  of  rotation 
B  =  sine  of  angle  of  rotation 
C  =  x,  y,  z  of  input  coordinates 
D  *  x,  y,  z  of  rotated  coordinates 


Subroutine  CSPR  (A,B,C) 
Program  CSPR 


This  subroutine  determines  the  Fresnel  integral  given  argument  C  by  call¬ 
ing  IBM  scientific  subroutine  CS. 

A  =  Real  part  of  Fresnel  integral  of  C 
B  ~  Imaginary  part  of  Fresnel  integral  of  C 
C  =  Argument  of  Fresnel  integral 


B.2  Ground  Reflection  Subroutines 


The  following  subroutines  related  to  the  ground  reflection  calculation 
(discussed  in  Section  3.2),  along  with  the  page  number  on  which  they  are  lo¬ 
cated,  are: 


Subroutine 

Page 

UtlKMLi 

j  s 

GREFC 

B-9 

VSUB 

B-10 

DTDEL 

B-10 

RCVANG 

B-10 

GRDOT 

B-i  1 

CRDRTB 

B-ll 

ROT 

B-l  1 

GRND 

B-12 

FREKFS 

B-l  2 

DOT 

B-l  3 

ODIST 

B-l  3 

DOTC 

B-l  3 

DREFC 

B-l  3 

SPECGD 

B-l  4 

CHSPEC 

B-l  3 

WGTAHP 

B-lb 

GXMTR 

B-l  7 

SETVAL 

B-l  7 

DIGLOD 

B-i  8 

GAIN 

B-19 

- t 

Subroutine  CHKMLT  (A) 

Program  HOJRDT 


See  Section  8.1  above. 
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Subroutine  GREFC 
Program  GREFC 


This  subroutine  determines  ground  reflection  parameters  if  IFGREF  is  1 
and  logical  variable  HUMP  is  false.  The  ground  plane  can  be  divided  into  10 
rectangular  elements  and  10  triangular  elements,  each  with  its  own  character¬ 
istic  roughness  coefficient  (SH2S)  and  dielectric  constant  (ERS).  The  remain¬ 
ing  regions  are  given  default  values  of  surface  roughness  (SH20)  and  dielec¬ 
tric  constant  (ERO).  Specular  reflection  parameters  for  each  of  these  regions 
if  ISPGRD  is  specified  as  0  in  block  data  BDUHFM.  If  ISPGRD  is  specified  as  1 
in  the  block  data,  a  simplified  method  is  used  to  compute  the  multipath  param¬ 
eters  as  noted  in  Section  3.2. 

INPUTS 

are  described  in  GREFC  comment  statements. 

OUTPUTS 

ARGUMENT  LIST  AMPL 

PHAS 
AZ 
EL 

TDEL 
AZI 
ELI 
OOP 
DOPR 


-  Relative  amplitude  of  ground  reflection 

■  Relative  phase  of  ground  reflection 

-  Planar  azimuth  in  inertial  XY  plane 

■  Planar  elevation  angle  in  inertial  XY 
plane 

■  Relative  time  delay  of  ground  reflection 

■  Planar  azimuth  receiver  incidence  angle 

■  Planar  elevation  receiver  incidence  angle 

■  Doppler  shift  of  ground  reflection 

-  Receiver  component  of  Doppler  shift 


Subroutine  VSUB  (A,B,C) 
Program  VSUB 


This  subroutine  determines  the  difference  between  two  vectors. 

A  •  An  array  of  three  components  for  Vector  i. 

B  -  An  array  of  three  components  for  Vector  2. 

C  -An  array  of  three  components  of  the  vector  dif¬ 

ference  between  vectors  1  and  2. 


Subroutine  DTDEL  (A,B,C,D) 
Program  DTDEL 


This  subroutine  determines  the  time  delay  of  a  specular  ground  reflection 
with  respect  to  the  direct  wave. 

A  “  Coordinates  of  the  3-component  vector  of  the 
transmitter  relative  to  specular  point  (in  iner¬ 
tial  space) 

B  “  Coordinates  of  the  3-component  vector  of  the  re¬ 
ceiver  relative  to  the  specular  point  (in  iner¬ 
tial  space) 

C  “  Coordinates  of  the  3-component  vector  of  the  re¬ 
ceiver  relative  to  the  transmitter  (in  inertial 
space) 

D  “  Time  delay  of  the  specular  ground  reflection 
relative  to  the  direct  wave 


Subroutine  RCVANG  (A,B,C,D,E,F,G,H) 
Program  GDDRONE 


See  Section  B.l  above. 
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This  subroutine  transforms  coordinates  from  the  inertial  system  to  that 
defined  by  the  Ath  surface  element. 

A  ■  Index  which  specifies  which  surface  element  is 
considered 

B  -  Original  coordinates  relative  to  Inertial  space  ; 

C  ■  Transformed  coordinates  relative  to  the  Ath  j 

surface  element  ! 

e 

l 


This  subroutine  transforms  coordinates  from  the  system  defined  by  the 
Ath  surface  element  of  the  inertial  system. 

A  -  Index  which  specifies  which  surface  element  is 
considered 

B  -  Coordinates  in  system  defined  by  the  Ath  sur¬ 
face  element 

C  ■  Coordinates  relative  to  inertial  coordinate  sys¬ 
tem 


Subroutine  ROT  (A,B,C,D) 
Program  ROT 


This  subroutine  applies  a  rotation  of  C  radians  about  the  D  axis  to  the 
vector  A. 

A  "  Original  coordinates 

B  -  Transformed  coordinates  assuming  rotation 
C  ■  Size  of  angle  (radians) 

D  ■  Indicates  which  axis 

B-ll 


Subroutine  GRDRTB  (A,B,C) 
Program  SPECGD 


Subroutine  GRDROT  (A,B,C) 
Program  SPECGD 


I 


1 


Subroutine  GRND  (A,B,C) 

Program  GRND 

This  subroutine  determines  which  surface  element  contains  a  given  point 
in  the  X-Y  plane. 

A  ■  Coordinates  of  integration  point 
B  -  Number  of  the  surface  element  which  contains  the 
point  under  consideration. 

C  *  Indicates  which  surface  elements  are  to  be 
searched 


Subroutine  FRERFS  (A,B) 
Program  FRERFS 


This  subroutine  computes  p8  which  is  the  product  of  the  Fresnel  reflec¬ 
tion  coefficient,  R,  and  the  attenuation  factor,  pr,  as  opposed  to  using  the 
Fresnel-Kirchof f  diffraction  formula  and  performing  a  detailed  confutation. 
This  subroutine,  as  noted  in  Section  3.2,  is  called  if  ISPGRD  is  specified  as 
1  in  block  data  file  BDUHFM. 

A  ■  Direction  vector  from  the  transmitter  to  the 
specular  point 

B  ■  Complex  specular  ground  reflection  coefficient 
Ps 

INPUTS 

LABELLED  COMMON 

XMTER  PLZ  -  True  — >  vertical  polarization 

-  False  — >  horizontal  polarization 

K  -  Wave  number 

K2  «  k2 
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LABELLED  COMMON 


SURF AC 

ERG  - 

Complex 

face 

SH2G  - 

Squared 

face 

dielectric  constant  of  ground  sur- 
RMS  roughness  height  of  ground  sur- 


Functlon  DOT  (A,B) 
Program  DOT 


This  function  computes  the  dot  product  of  two  vectors  A  and  B. 


Function  DDIST  (A) 
Program  DDIST 


This  function  computes  the  length  of  a  vector  with  coordinates  A  using 
double  precision. 


Function  DOTC  (A,B,C,D,E,F) 
Program  DOTC 


See  Section  B.l  above. 


Function  DREFC  (A,B) 
Program  DREFC 


This  subroutine  calculates  the  differential  reflectivity  coefficient  for 
a  given  grid  point  used  in  the  approximation  of  the  Fresnel  integral. 
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A  ■  3-cooponent  vector  giving  the  coordinates  of  the 
surface  eleaent 

B  -  Index  of  surface  eleaent  containing  point  A. 


Subroutine  SPECGD  (A,B,C,D) 
Program  SPECGD 


This  subroutine  confutes  the  location  of  the  specular  point  for  an  arbi¬ 


trary  tilted  ground  region. 

A 

B 

C 

D 


Location  of  specular  point 


X 


on 


In 


Parameters  specifying  the  size 
size  of  the  region  over  which 
the  Fresnel  integration  is  to 
be  performed. 


KEY  INPUTS 
LABELLED  COMMON 

GRDRT  RECTA  ■  X,Y  coordinates  of  fourth  corner  of  I-th 

rectangle  (found  in  subroutine  INIT) 

RTCRS  ■  Orientation  vector  for  I-th  rectangle  or 
triangle  (found  in  subroutine  INIT) 


RTSPEC 


GRDTRN 


MRTSPC  •  Indices  to  be  used  in  computing  specular 
ground  reflection  for  rectangles  and  trian¬ 
gles  (found  in  subroutine  INIT) 

NSPEC  -  Number  of  specular  ground  reflection  compo¬ 
nents  (found  in  subroutine  INIT) 


DICS  ■  9  direction  cosines  to  specify  the  orienta¬ 

tion  of  the  surface  elements 
XYZOG  ■  Coordinates  to  translate  origin  for  each 
surface  eleaent  considered 
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SLOPE 

- 

Slopes  of  sides  of  triangles 

YINT 

- 

Y  Intercept  of  sides  of  each  triangle 

XMTER 

XMTR 

am 

Transmitter  position 

XVEL 

- 

Transmitter  velocity 

RCVER 

RCVR 

m 

Receiver  position 

SURFACE 

NS 

- 

Number  of  surface  elements  (less  than  or 

equal  to  20) 

Subroutine  CHSPEC  (A,B,C) 

Program  SPEOGD 

This  subroutine  repositions  the  specular  point  for  the  surface  ele¬ 
ment 

A  •  Index  to  specify  which  surface  element  is  being 
considered. 

B  ■  Coordinates  of  specular  point. 

C  *  Coordinates  of  repositioned  specular  point. 


INPUTS 

LABELLED  COMMON 

SURF AC  NT  «  Number  of  triangles 

NR  -  Number  of  rectangles 

SURFS  ■  Defining  coordinates  of  the  surface  ele 
meats 


B-15 


Subroutine  UCTAMP  (A,B,C,D,E) 
Program  MLTPLOT 


This  subroutine  applies  amplitude  weighting  due  to  the  transmitter  an¬ 
tenna  pattern. 

A  *  Multiple  emitter  index 

B  -  Adjusted  multipath  amplitude  normalized  with 
respect  to  that  of  the  direct  wave  from  the  azi¬ 
muth  transmitter 

C  -  Multipath  azimuth  angle  at  the  transmitter 
D  <*  Multipath  elevation  angle  at  the  transmitter 
E  *  Raw  multipath  amplitude 


KEY  INPUTS 
LABELLED  COMMON 


AZBLNK 

NA2BLK, 

tazblk 

DMBLNK 

NDMBLK, 

TDMBLK 

EIBLNK 

NE1BLK, 

TE1BLK 

E2BLNK 

NE2BLK, 

NE2BLK 

where 

N _ BLK 

“ 

Number  of  off-time  intervals 

T _ BLK 

- 

Off-time  intervals 

DIRWAV 

AZD1R 

m 

Planar  azimuth  angle  of  direct  wave 

ELDIR 

m 

Planar  elevation  angle  of  direct  wave 

CMFLTP 

TIME 

- 

Time  along  flight  path 

BXMPAR 

KXMTRK 

m 

Tracking  option  parameter 

BAZTRK 

m 

Planar  azimuth  pointing  angle 

BELTRK 

- 

Planar  elevation  pointing  angle 
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Function  GXMTR  (A,B,C) 

Program  JCHOJ 

This  subroutine  is  used  to  determine  the  gain  of  the  A**1  transmitter 
given  azimuth  and  elevation  planar  angles  of  B  and  C.  This  subroutine  is 
called  twice.  The  first  time  modifications  to  the  specular  reflection  param¬ 
eters  are  determined.  The  second  time  adjustments  are  made  to  the  direct 
wave. 


Subroutine  SETVAL  (A,B,C,D,E,F,G,H,1,J 
Program  SETVAL 


This  subroutine  is  called  to  transfer  parameters  associated  with  direct 
wave  and  specular  reflection  to  labelled  common  areas  MULTI  and  MULTIP. 

A  •  Index  indicating  where  quantities  are  to  be 
loaded  in  labelled  common  arrays. 

B  *  Amplitude 
C  “  Phase 

0  “  Azimuth  planar  angle  relative  to  the  transmitter 

E  ■  Elevation  planar  angle  relative  to  the  transmit¬ 
ter 

F  ■  Time  delay 

G  *  Magnitude  of  the  Doppler  shift 
H  “  Receiver  component  of  the  Doppler  shift 
I  "  Azimuth  planar  angle  relative  to  receiver 
J  ■  Elevation  planar  angle  relative  to  receiver 
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Subroutine  DIGLOD  (A,B,C,D,E,F,G,H,I,J) 
Program  MLTPLOT 


This  subroutine  Is  called  to  transfer  parameters  associated  with  the  di¬ 
rect  wave  and  specular  reflection  to  the  labelled  common  area  DIAGP. 

A  -  Index  Indicating  which  quantitites  are  to  be 
used  in  labelled  common  MULTI. 

B  -  Signal  processing  Index  (2  for  the  enhanced  In¬ 
terferometer) 

C  -  Multiple  emitter  index 

D  •  Index  specifying  which  multipath  obstacle  is  un¬ 
der  consideration  (1  for  the  direct  wave,  2  for 
the  ground,  3  through  12  for  buildings,  13 
through  22  for  aircraft). 

E  ”  Weighted  multipath  amplitude  normalized  with 
respect  to  that  of  the  direct  wave  from  emitter 
1. 

F  »  Raw  amplitude 

G  -  Multipath  planar  azimuth  angle  relative  to  emit¬ 
ter  axes 

H  «  Multipath  planar  elevation  angle  relative  to 
emitter  axes 

I  ■  Receiver  fractional  Doppler  frequency 

J  ■  Time  delay 
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Subroutine  GAIN  (A,B,C,D,E) 

Program  CCGAIND 

This  subroutine  characterizes  the  receiver  antenna  pattern: 

A  -  Index  of  multipath  component 
B  ■  Antenna  element 
C  *  Emitter  number 
D  -  Modified  phase 
E  -  Modified  gain 

B.3  Building  Reflection  Subroutines 

The  following  subroutines  related  to  building  reflectior 
(discussed  in  section  3.3),  along  with  the  page  number  on  whict 
cated,  are: 


Subroutine 

CHKMLT 

BREFC 

B COORD 

BCOMP 

CSPR 

VSUB 

DTDEL 

BCOORT 

POLZN 

RCVANG 

VNORM 

VNORMC 

WGTAMP 

GXMTR 

SETVAL 

DIGLOD 

GAIN 


Page  Number 
B-20 
B-20 
B-20 
B-21 
B-21 
B-21 
B-21 
B-22 
B-22 
B-23 
B-23 
B-23 
B-23 
B-23 
B-24 
B-24 
B-24 


Iculatlons 
are  lo- 


Subroutine  CHKMLT  (A) 
Progr sp  ’':JRDT 


See  Section  B.l  above. 


Subroutine  BREFC  (A, B ,C ,D ,E ,F ,G ,H,I , J ,K) 
Program  BREFC 


This  subroutine  determines  the  reflection  multipath  for  a  given  building. 

A  *  Index  specifying  which  building  is  being  consid¬ 
ered 

B  »  Index  specifying  which  path  is  being  considered 
C  »  Relative  amplitude  of  the  building  reflection 
D  =  Phase  of  the  building  reflection 
E  =*  Planar  azimuth  angle 
F  =  Planar  elevation  angle 

G  =  Relative  time  delay  of  the  building  reflection 
H  =  Doppler  shift  of  the  building  reflection 
I  =  Receiver  component  of  building  reflection  Dopp¬ 
ler  shift 

J  »  Planar  incident  azimuth  angle 
K  «*  Planar  incident  elevation  angle 


Subroutine  BCOORD  (A,B,C) 

Program  BCOORD 

This  subroutine  converts  ground  position  coordinates  to  a  building  cen¬ 
tered  coordinate  system. 

A  -  Building  number 

B  =  Ground  coordinates  of  positioi,  vector 
C  -  Building  centered  coordinates  of  position  vector 
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Entry  BCOMP  (A,B,C) 
Program  BCOORD 


This  subroutine  converts  velocity  vectors  in  a  ground  system  to  a  build¬ 
ing  referenced  system. 

A  -  Building  number 

B  *  Ground  coordinates  of  velocity  vector 
C  -  Building  centered  coordinates  of  velocity  vector 


Subroutine  CSPR  (A,B,C) 
Program  CSPR 


Fresnel  Integral  subroutine.  See  Section  B.l. 


Subroutine  VSUB  (A,B,C) 
Program  VSUB 


Vector  subtraction  subroutine.  See  Section  B.2 


Subroutine  DTDEL  (A,B,C,D) 
Program  DTDEL 


Subroutine  for  determining  time  delay  associated  with  a  building  multi- 
path  relative  to  the  direct  wave.  See  Section  B.2 
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Subroutine  BCOORT  (A,B,C) 

Program  BCOORT 

This  subroutine  converts  building  referenced  coordinates  of  a  position 
vector  to  a  ground  referenced  coordina  a  system. 

A  -  Building  number 

B  -  Building  coordinates  of  position  vector 
C  -  Ground  coordinates  of  position  vector 


Subroutine  POLZN  (A,B,C,D,E,F,G) 
Program  POLZN 


This  subroutine  computes  the  equivalent  Fresnel  reflection  coefficient 
for  reflection  from  a  tilted  planar  facet. 

A  -  Relative  complex  dielectric  constant 
B  *  Ground  unit  normal  vector  (not  used  in  calcula¬ 
tions) 

C  -  Planar  facet  unit  normal  vector 
D  -  Vector  from  specular  point  to  the  transmitter 
E  -  Vector  from  specular  point  to  the  receiver 
F  ■1  for  a  perfect  conductor;  0  otherwise 
G  «  Equivalent  Fresnel  reflection  coefficient 

KEY  INPUT 
LABELLED  COMMON 

XMTR  PLZ  -  Logical  variable  describing  emitter  polari¬ 

zation;  true  for  vertical  polarization, 
false  for  horizontal  polarization. 
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See  Section  B.2 


Subroutine  SETVAL  (A,B,C,D,E,F,G,H,I, J) 
Program  SETVAL 


See  Section  B.2. 


Subroutine  DIGLOD  (A,B,C,D,E,F,G,U,I, J) 
Program  MLTPLOT 


See  Section  B.2. 


Subroutine  GAIN  (A,B,C,D,E) 
Program  CCGAIND 


See  Section  B.2. 


B.4  Aircraft  Reflection  Subroutines 

The  following  subroutines  related  to  building  reflection  calculations 
(discussed  in  Section  3.4),  along  with  the  page  number  on  which  they  are 
found,  are: 


Subroutine 

CHKMLT 

FREFC 

ACTYPE 

ORIENT 

CYLREF 

VSUB 

XDIV 


Page  No. 
B-25 
B-25 
B-26 
B-27 
B-27 
B-28 
B-28 


B-24 


r 


NORMLZ 

B-28 

ROTATE 

B-29 

CSPR 

B-29 

ORIEND 

B-29 

POLZN 

B-30 

TREFC 

B-30 

WGTAMP 

B-30 

GXMTR 

B-30 

SETVAL 

B-30 

D1GL0D 

B-3 1 

GAIN 

B-31 

MAXMLT 

B-3 1 

PTIME 

B-31 

DIGOUT 

B-31 

Subroutine  CHKMLT  (A) 
Program  HOJRDT 


See  Section  B.l  above. 


Subroutine  FREFC  (A,B,C,D,E,F,G,H,I,J,K) 
Program  FREFC 


This  subroutine  determines  the  multipath  quantities  associated  with  re¬ 
flections  off  aircraft  fuselages. 

A  •  Airplane  number 

B  ■  Multipath  component  number 

C  -  Relative  amplitude  of  fuselage  reflection 
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D  ■  Relative  phase  of  fuselage  reflection 
E  ■  Planar  azimuth  angle  with  respect  to  the  trans¬ 
mitter 

F  -  Planar  elevation  angle  with  respect  to  the 
transmitter 

G  ■  Relative  time  delay  of  fuselage  reflection 
H  ■  Doppler  shift  of  fuselage  reflection 
I  ■  Receiver  component  of  Doppler  shift 
J  »  Planar  azimuth  angle  with  respect  to  receiver 
K  ■  Planar  elevation  angle  with  respect  to  receiver 


Subroutine  ACTYPE  (A,B,C,D,E,F,G,H) 
Program  JCHOJ 


This  subroutine  sets  the  parameters  for  specified  aircraft  which  charac¬ 
terizes  their  geometry.  These  quantities  are  presently  specified  for  aircraft 
such  as  the  DC-10  and  the  707  because  the  program  was  originally  designed  to 
study  multipath  environments  in  the  vicinity  of  airfields. 

A  »  Aircraft  type 

B  ■  Fuselage  radius 

C  »  Fuselage  length 

D  ■  Tail  fin  radius  of  curvature 

E  ■  Length  of  tail  fin 

F  *  Height  of  tail  fin 

G  *  Total  area  of  both  wings 

H  ■  Height  of  center  of  fuselage  above  ground 
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Subroutine  ORIENT  (A,B,C) 
Program  ORIENT 


This  subroutine  confutes  the  equivalent  vertical  cylindrical  components 
of  a  vector  given  Its  ground  components  with  respect  to  a  horizontal  cylindri¬ 
cal  fuselage. 

A  -  Number  specifying  which  airplane  fuselage  Is  un¬ 
der  consideration. 

B  -  Ground  components  of  vector  with  respect  to  a 
horizontal  cylinder 
C  »  Vertical  cylindrical  components 


Subroutine  CYLREF  (A,B,C,D,E,F,G,H,I,J,K,L) 
Program  CYLREF 


This  subroutine  computes  the  multipath  parameters  of  a  smooth,  perfectly- 
conducting  section  of  a  vertical  cylindrical  surface. 

A  -  Index  to  specify  which  airplane  is  under  consid¬ 
eration 

B  »  Index  to  specify  which  multipath  component  is 
under  consideration 

C  *  Reflection  amplitude 

D  ■  Reflection  phase 

E  -  Planar  azimuth  angle  with  respect  to  inertial  X 
axis 

F  -  Planar  elevation  angle  with  respect  to  Inertial 
X  axis 

G  -  Relative  time  delay 

H  -  Doppler  shift 
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I  -  Receiver  component  of  Doppler  shift 
J  ■  Specifies  which  vertical  direction  represents 
true  vertical  (so  as  to  use  the  same  program  for 
modelling  reflections  from  horizontal  fuselages 
and  vertical  tail  fins) 

K  ■  Planar  incident  azimuth  angle 
L  -  Planar  incident  elevation  angle 


Subroutine  VSUB  (A,B,C) 
Program  VSUB 


See  Section  B.2. 


Subroutine  XDIV  (A,B) 
Program  NORMLZ 


This  is  the  entry  point  to  subroutine  NORMLZ  used  to  normalize  a  vector. 

A  -  Coordinates  of  vector 
B  -  Length  of  vector 


Subroutine  NORMLZ  (A) 
Program  NORMLZ 


This  subroutine  normalizes  a  vector. 

A  -  Coordinates  of  vector. 


B-28 


Subroutine  ROTATE  (A,B,C) 
Program  ROTATE 


This  subroutine  finds  the  conponents  of  a  dimensional  vector  in  a  rotated 
coordinate  system. 

A  -  Coordinates  of  vector  in  unrotated  coordinate 
system 

B  -  Coordinates  of  vector  in  rotated  coordinate  sys¬ 
tem 

C  -  Coordinates  of  unit  vector  in  rotated  x  direc¬ 
tion  relative  to  unrotated  system. 


Subroutine  CSPR  (A,B,C) 
Program  CSPR 


See  Section  B.l. 


Subroutine  ORIEND  (AtB,C) 
Program  P0U5N 


This  subroutine  transforms  the  vector  coordinates  from  the  cylinder  coor¬ 
dinate  system  to  the  lnertially  referenced  system. 

A  -  Airplane  number 

B  -  Vector  in  fuselage  coordinate  system 
C  -  Vector  in  ground  coordinate  system 
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Subroutine  POLZN  (A,B,C,D,E,F,G) 
Program  POLZN 


See  Section  B.3. 


See  Section  B.3 


Subroutine  GAIN 


This  subroutine  determines  the  scaling  parameters  for  multipath  diagnos¬ 
tics  (plots)  at  a  particular  evaluation  point.  The  quantities  determined  are 
Loaded  into  labelled  common  DIGSCL. 


See  Section  B.3. 


Subroutine  MAXMLT  (A) 
Program  MAXMLT 


Subroutine  PTIME  (A,B,C,D,E,F,G) 
Program  PTIME 


This  subroutine  is  used  to  determine  the  CPU  time  required  for  multipath 
calculations.  See  Section  2.0  and  Appendix  A. 


Subroutine  OIGOUT  (A,B) 
Program  MLTPLOT 


This  subroutine  writes  out  multipath  data  to  disk  (UNIT  14) 


Appendix  C 

SUBROUTINES  FOR  SEEKER  MODULE 


This  appendix  provides  additional  details  of  the  subroutines  used  for  the 
seeker  calculations.  It  is  most  useful  when  it  is  used  with  the  listings 
Included  with  this  report. 

The  following  subroutines  are  Included  along  with  the  page  number  of 
which  they  are  found. 


Subroutine 

Page  No 

HOJ 

C-3 

RDPHAS 

C-3 

DLS 

C-3 

Cl  HOJ 

C-4 

GAIN 

C-4 

RESOLV 

C-5 

VSUB 

C-5 

CONVER 

C-5 

ERRMOD 

C-5 

EMTLOD 

C-6 

UPDTMS 

C-6 

ATRPLT 

C-6 

ATRPL2 

C-6 

TRACK 

C-7 

EMCASE 

C-7 

RCVANG 

C-7 

INVERT 

C-8 

ATRPLT 

C-8 

PUTHOJ 
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C.l  GENERAL  REQUIREMENTS  AND  DEFINITIONS 

Requirements  for  Interferometer  Signal  Processing  to  be  Used 

1.  Specify  TDG  in  response  to  question  on  signal  processing  op¬ 
tion  to  be  used  in  exec  TESTSIM2. 

2.  Load  CAT1,  CSRPVTRK,  C1H0JG2,  and  CCGAIND  in  that  order  in 
the  load  sequence 
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Quantities  Required  for  Signal  Processing  Computations 

The  following  quantities  are  to  be  specified  by  the  user  in  the  block 
data  file  BDUHFM: 


FAD2  - 
FED2  - 
F1A  - 
FIE  - 


DYNX  * 

AMI  5  » 

ITRK  - 
IGAIN  - 
BUT  - 

ORATE  - 


NXMTR  - 


(Azimuth)  d/2 
(Elevation)  d/2 
(Azimuth)  2ifd/X 
(Elevation  2*d /X 


(Meters) 

(Meters 

(Dimensionless) 

(Dimensionless) 


where  d  is  the  spacing  between  antenna  elements,  X  is 
the  wavelength  of  the  incident  radiation 
Scaling  parameter  X  in  the  relationship 
9  -  X(PHIP)  (Default  -  1) 

Quantity  used  for  resolving  phase  difference  (Default  ■ 

0) 


1  for  tracker  implementation 

1  for  modelling  of  seeker  antenna  characteristics 
Coefficient  used  in  some  version  of  the  subroutine  mod¬ 
elling  antenna  gain  (although  not  in  CCGAIND). 
Simulation  interval  for  determining  signal  in  space  pa¬ 
rameters,  determining  planar  angle  estimates,  and  up¬ 
dating  autopilot  (nominally  0.1  sec.  for  this  vehicle) 
Number  of  emitters  (less  than  or  equal  to  4) 


The  following  quantities  are  determined  in  the  multipath  computation 
b lock : 


AMP  - 
PHASE  - 
TAU  - 
RDOP  - 
AZP  - 


multipath  amplitude 
multipath  phase 
multipath  time  delay 

multipath  fractional  doppler  frequencies 
multipath  azimuth  incidence  angles  relative  to  the  re¬ 
ceiver 
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ELP  “  multipath  elevation  incidence  angles  relative  to  the 
receiver 

OMEG  -  frequency  of  emitter  signal 

C .2  SEEKER  MODULE  SUBROUTINES 


Subroutine  HOJ  (IANT,  MISSLE) 
Program  HOJRDT 


This  subroutine  is  the  driver  which  determines  the  planar  angles  that 
specify  the  orientation  of  the  estimated  line-of-sight  vector  with  respect  to 
the  vehicle  boresight.  When  the  enhanced  interferometer  is  to  be  invoked  then 
the  parameter  MISSILE  is  2  and  the  parameter  IANTI  is  set  to  1  in  the  overall 
simulation  driver,  JCMLST. 


j 

Subroutine  RDPHAS  (IY) 

Program  HOJRDT 

Specifying  TDG  in  the  EXEC  file  TESTSIM2  as  noted  in  A.l  calls  this  sub¬ 
routine  which  adds  a  random  phase  angle  between  0  and  2ir  to  the  multipath 
phase  variable  phase  contained  in  labelled  common  MULTI.  This  subroutine 
calls  a  system  subroutine  for  generating  the  random  number  IY. 


Subroutine  DLS  (ITH) 
Program  CAT1 


As  for  RDPHAS,  specifying  TDG  in  the  EXEC  calls  this  subroutine  which 
calls,  in  turn,  the  subroutine  CIHOJ  for  the  enhanced  interferometer.  This 
subroutine  generates  the  planar  angles  by  modelling  the  antenna  and  signal 
processor  end  the  tracker  if  it  is  invoked. 
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Note  that  RDPHAS  and  DLS  are  called  twice.  On  the  first  call,  the  azi¬ 
muth  planar  angle  Is  calculated  and  on  the  second,  the  elevation  angle.  The 
following  array  contains  the  angles: 

HOJERR(l)  «  Azimith  planar  angle 

H0JERR(2)  -  Not  used 

H0JERR(3)  -  Elevation  planar  angle 

This  subroutine  also  generates  planar  angles  between  the  vehicle  velocity  vec¬ 
tor,  boresight  vector,  and  seeker  estimated  LOS  relative  to  the  true  line-of- 
sight  to  the  emitters  in  the  scenario  that  has  been  specified. 


Subroutine  C1H0J  (IP) 
Program  CLH0JG2 


This  subroutine  models  the  antenna  pattern,  signal  and  angle  processing. 
The  parameter  IP  is  set  to  1  in  the  calling  subroutine  DLS.  Details  of  this 
subroutine  are  given  in  Section  4.2. 


Subroutine  GAIN  (A,B,C,D,E) 
Program  CCGAIND 


This  subroutine  models  the  specific  characteristics  of  the  EHT  antenna 
given  the  angle  of  arrival  of  the  multipath  signals.  The  parameters  passed 
are: 

A  *  Index  of  the  multipath  signal 
B  =  Antenna  element 
C  *  Index  of  the  emitter 
D  *  Phase  change 
E  =*  Antenna  gain 
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Subroutine  RESOLV  (A,B,C,D) 
Program  C1H0JG2 


This  subroutine  resolves  angle  A  into  the  angular  region  defined  by  an¬ 
gles  B  and  C.  The  output  of  the  subroutine  is  angle  D. 


Subroutine  VSUB  (A,B,C) 
Program  VSUB 


This  subroutine  computes  the  vector  difference  between  A  and  B  giving  the 
result  in  C. 


Subroutine  CONVER  (A,B,C) 
Program  CIKOJG2 


This  subroutine  converts  the  planar  azimuth  and  elevation  angles  associ¬ 
ated  with  the  Ath  multipath  component  to  the  conical  angles  B  and  C. 


Subroutine  ERRMQD  (A,B) 
Program  ERRMOD 


This  subroutine  reads  input  defining  the  magnitude  and  characteristics  of 
error  sources  and  perturbs  the  angles  determined  by  the  enhanced  interferome¬ 
ter  processor.  Gaussian  noise  corrupting  the  angles  is  modelled  as  well  as 
constant  biases  and  blind  range  effects  as  well  as  errors  associated  with  the 
deviation  from  the  ideal  interferometer  relationship.  These  errors  are  based 
on  measurements.  The  parameters  are  defined  as  follows: 
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A  “  Perturbed  azimuth  angle  (radians) 

B  -  Perturbed  elevation  angle  (radians) 


Subroutine  EKTLOD  (A,B,C,D) 
Program  CRSPLT1 


This  subroutine  loads  the  x,y,z  coordinates  of  the  Atl*  emitter  Into  ar¬ 
guments  B,  C,  and  D. 


Subroutine  UPDTMS  (A,B,C) 
Program  UPDTMS 


This  subroutine  determines  the  miss  distance,  A,  from  the  Bctl  emitter 
given  the  vector  C  from  the  vehicle  to  that  emitter. 


Subroutine  ATRPLT  (A,B,C,D,E) 
Program  GDDRONE 


This  subroutine  interpolates  tables  A  and  B  to  find  E  for  a  given  D. 
These  tables  have  C  element. 


Subroutine  ATRPL2  (A,B,C,D,E,F,G,H) 
Program  GDDRONE 


This  subroutine  performs  a  two  dimensional  interpolation  using  tables  A 
and  B  with  dimensions  D  and  E,  respectively,  by  looking  up  values  in  table  C. 
Arguments  F  and  G  are  used  for  determination  of  the  proper  indices.  The  out¬ 
put  is  H. 
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Subroutine  TRACK  (A,B) 
Program  CSRPVTRK 


This  subroutine,  which  Is  described  In  some  detail  In  Section  A. 3,  takes 
planar  azimuth  and  elevation  angles  A  and  B  and  replaces  them  If  necessary 
with  revised  values  based  upon  several  criteria. 


Subroutine  EMCASE  (A,B) 
Program  CRSPLT1 


This  subroutine  determines  the  vector  B  from  the  vehicle  to  the  Atl1 
emitter. 


Subroutine  RCVANG  (A,B,C,D,E,F,G,H, ) 
Program  GDDRONE 


This  subroutine  determines  azimuth  and  elevation  planar  angles  In  the  re¬ 
ceiver  coordinate  system  given  the  inertially  referenced  coordinates  of  the 
vector  from  receiver  to  target. 

A  -  X  ) 

vector  coordinates  to  a  point  in 

B  -  Y  } 

I  space  in  inertial  reference  frame 

C  -  Z 
D  =  dummy 
E  »  dummy 
F  ■  dummy 

G  *■  planar  azimuth  angle  relative  to  vehicle  boresight 
H  -  planar  elevation  angle  relative  to  vehicle  boresight 
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Subroutine  INVERT  (A,B,C) 
Program  GDDRONE 


Thia  subroutine  determines  the  vector  coordinates,  C,  in  t  e  ne 
frame  using  vector  coordinates,  B,  in  the  body  frame  and  coordinates  o 
boresight  vector  A  in  the  inertial  frame. 


Subroutine  ATRPLT  (A.B.C.D.E) 

Program  ERRMOD 

This  subroutine  interpolates  in  tables  A  and  B  to 
The  tables  have  C  elements. 


find  E  for  a  given  C. 


Subroutine  PUTHOJ  (A,B) 

Program  CRSPLT2 

This  subroutine  writes  out  the  vehicle  position, 
coordinates,  planar  angle  estimates  from  the  tracker, 
difference  angles  to  the  disk. 


velocity  and  boresight 
error  angles  and  planar 
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Appendix  D 

SUBROUTINES  FOR  FLIGHT  DYNAMICS  MODULE 


This  appendix  provides  additional  details  relative  to  the  subroutines 
used  by  the  flight  dynamics  module.  It  is  most  useful  when  used  with  the 
listings  included  with  this  report. 

The  following  are  the  subroutines  included  along  with  the  page  number  on 
which  they  are  found. 


FLTDYN 

CRSGID 

AUTPLT 

AERDYN 

EQAMO 

KINEM 

WIND 

RCWIND 

RCVANG 

MICRO 

ATRPL2 

BLIMIT 

INTGRT 

INVERT 

GYERR 


Page  No 

D-l 
D-2 
D-3 
D-4 
D-6 
D-8 
D-9 
D-IO 
D-10 
D-l  1 
D-l  1 
D-l  2 
D-12 
D-l  3 
D-13 


Subroutine  FLTDYN  (A) 
Program  GDDRONE 


This  subroutine  is  the  driver  which  controls  the  sequencing  of  all  of  the 
subroutines  in  the  flight  dynamics  module.  When  the  enhanced  Interferometer 
is  used  (which  is  the  case  for  the  description  given  in  this  report),  the 
parameter  A  is  set  to  2. 
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i 


INPUTS 

LABELLED  COMMON 


FLTIND 

INDEX 

Integer  Index  which  specifies  how 

many  interations  have  been  made 

in  the  multipath  computations. 

AIRPRT 

DRATE 

m 

Interval  over  which  multipath 

quantities  are  determined. 

OUTPUTS 

LABELLED  COMMON 

SAMHLD 

INTRTE 

m 

Interval  for  performing 

integration  of  the  equations  of 

motion. 

TWD 

m 

^  1 .  constants  associated  with  the 

WAV 

m 

*  1.  tunable  integration  scheme 

This  subroutine  Is  called  by  the  master  driver  for  the  simulation  SJCMLST. 


Subroutine  CRSGID  (PTCHCM,  YAWCM) 
Program  GDGUID 


This  Is  the  guidance  computer  subroutine  which  processes  the  planar  azi¬ 
muth  and  elevation  angles  to  produce  pitch  and  yaw  commands  for  the  auto¬ 
pilot.  The  inputs  and  outputs  of  the  subroutine  are: 


INPUTS 

LABELLED  COMMON 

H0J3  HOJERR(l)  -  Planar  azimuth  angle  (rads) 

H0JERR(3)  ■  Planar  elevation  angle 

(rads) 
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f 


GYRFD 

OAQ 

Filtered  value  of  corrupted  pitch 

and  body  rate  (deg/sec) 

AOR 

m 

Filtered  value  of  corrupted  yaw 

body  rate  (deg/sec) 

OUTPUTS 

SUBROUTINE 

ARGUMENTS 

PTCHCM 

m 

Pitch  steering  command  (deg/sec) 

YAWCM 

* 

Yaw  steering  command  (deg/sec) 

This 

subroutine 

is  called 

by  FLTDYN  the  driver  program  for  the  flight  dynamics 

part 

of  the  simulation. 

It  is  called 

every 

DRATE  seconds. 

Subroutine  AUTPLT  (PTCHCM,  YAWCM,  DELEN,  DELRN,  AILDEG) 

Program  GDDRONE 

This  is  the  autopilot  subroutine  which  accepts  commands  from  the  guidance 
computer  to  produce  elevator,  rudder,  and  aileron  positions.  The  inputs  and 
outputs  are: 

INPUTS 


SUBROUTINE  ARGUMENTS 

PTCHCM 

* 

Guidance  pitch  command 

YAWCM 

m 

Guidance  yaw  command 

LABELLED  COMMON 

GYRFD 

0AQ1 

m 

Corrupted  pitch  rate  measurement 

(after  modelling  of  the  gyroscope 

errors) 

0AR1 

Corrupted  yaw  rate  measurement 

(after  modelling  of  the  gyroscope 

errors) 

D-3 


H0J2 


P0SITN(3) 


Perfect  elevation  measurement 


DATRNS 

PSIN 

* 

Perfect  yaw  measurement 

XANROT 

FIRF 

- 

Perfect  roll  measurement 

OUTPUTS 

SUBROUTINE  ARGUMENTS 

DELEN 

- 

Elevator  deflection  (degs) 

DELRN 

- 

Rudder  deflection  (degs) 

AILDEG 

- 

Aileron  deflection  (degs) 

This  subroutine  is  called  by  FLTDYN  the  driver  program  for  the  flight  dynamics 
part  of  the  simulation  every  1NTRTE  seconds. 


Subroutine  AERDYN  (RALPHA,  RBETA,  DELEN,  DELRN,  AILDEG,  ALAPRAT, 
THREXT,  CXN,  CYN,  CZN,  COEFL,  COEFH,  COEFN,  VWXN,  VWYN,  VWZN,  SMP, 
SMQ,  SMR) 

Program  GDDRONE 


This  is  the  aerodynamics  subroutine  which  provides  the  aerodynamic  forces 
for  the  equations-of-motlon.  The  inputs  and  outputs  are: 


INPUTS 

SUBROUTINE  ARGUMENTS 


RALPHA 

m 

Angle  of  attack,  (rads) 

RBETA 

m 

Sideslip  (rads) 

DELEN 

m 

6E  (degs) 

DELRN 

m 

6r  (degs) 

AILDEG 

m 

6a  (degs) 

ALPRAT 

m 

Rate  of  change  of  a  (rads/sec) 
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VWXN 

VWYN 

VWZN 


Linear  and  angular  velocltie8 
associated  with  the  wind  field 


LABELLED  COMMON 
SWMARK 


SMP 

SMQ 

SMR 


PMODl  -  Logical  variable  to  indicate 
pitchover  (.TRUE.  implies  the 
vehicle  is  in  the  dive) 


DATRNS 

DELRNP 

tm 

SR 

(rads) 

DELENP 

m 

SE 

(rads) 

hiOJ2 

P0SITNO) 

m 

Vehicle  altitude 

AFRDBK 

OAQN 

m 

Q  (deg/sec) 

OARN 

m 

R  (deg/sec) 

OAPN 

m 

P 

(deg/sec) 

OUTPUTS 


CXBN 

m 

Cx 

CYN 

m 

Cy 

CZN 

m 

cz 

Aerodynamic 

COEFL 

m 

cl 

Coefficients 

COEFM 

m 

CM 

COEFN 

a* 

CN 

/ 

THREXT 

* 

Thrust 

This  subroutine  is  called  by  FLTDYN  the  driver  program  for  the  flight  dynamics 
part  of  the  simulation  every  1NTRTE  seconds. 
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Subroutine  EQUAMO  (THREXT,  CXN,  CYN,  CZN,  COEFL,  COEFN,  VELP 
RALPHA,  RBETA,  ALPRAT ,  VWXN,  VWYN,  VWZN,  SMP,  SMQ,  SMR) 
Program  GDDRONE 


I _ 

This  Is  the  equatlons-of -motion  subroutine  which  computes  the  accelera¬ 
tions  and  velocities  of  the  vehicle  in  its  own  axis  system,  along  with  the 
space  orientation.  In  addition,  the  velocities  of  the  wind  field  relative  to 
the  body  axes  and  corruption  of  the  body  angular  rates  by  gyro  errors  are  also 
computed.  The  inputs  and  outputs  are: 


INPUTS 

SUBROUTINE  ARGUMENTS 


OUTPUTS 

LABELLED  COMMON 
ARFDBK 


THREXT  -  Thrust 

CXN  -  Cx 

CYN  -  Cy 

CZN  -  C2 

COEFL  *  CL 

COEFM  -  CM 

COEFN  -  CN 


OAQN  ■  Q  (degs/sec) 
OARN  “  R  (degs/sec) 
OAPN  -  P  (degs/sec) 


VUN 

-  vx; 

velocity 

along 

the 

body 

X 

axis 

(ft/sec) 

VVN 

-  VY; 

velocity 

along 

the 

body 

y 

axis 

(ft/sec) 

VWN 

“  VZ; 

velocity 

along 

the 

body 

z 

axis  (ft/sec) 


D-6 


ARGUMENT 


LABELLED 

DATRNS 


H0J2 


ROLPLT 


XANROT 


LIST 


COMMON 


VELP(3) 

- 

VELP(l)  - 

VELP(2 )  - 

Vy 

VELP(3)  - 

VZ 

RALPHA 

4 

m 

a  (rads/sec) 

RBETA 

m 

8  (rads) 

ALPRAT 

m 

• 

a  (rads/sec) 

V 

linear  and 

VWXN 

m 

vwx  (ft/sec) 

angular 

VWYN 

m 

Vwy  (ft/sec) 

velocities 

VWZN 

m 

VW2  (ft/sec) 

> 

determined  in 

SMP 

m 

Pw  (rads/sec) 

the  call  to 

SMQ 

as 

Qy  (rads/sec) 

subroutine 

SMR 

m 

Ry,  (rads/sec) 

UIND 

ATACKP 

m 

a  (degs) 

NPN 

m 

Pitch  load  factor  (ft/sec2) 

ATACKY 

m 

6  (degs) 

NYN 

m 

Yaw  load  factor  (ft/sec2) 

AQNP 

as 

Q  (rads/sec) 

THETN 

m 

0  (rads) 

ARNP 

as 

R  (rads/sec) 

PSIN 

■s 

ip  (rads) 

BOREST(3) 

m 

Projection  of  unit 

pointed  along  vehicle  bore- 
sight  on  inertial  X,Y,Z  axes 


FI  DOT 

m 

* 

(rads/sec) 

FI 

- 

* 

(rads  where 

-0°  <  $  <  ®) 

FIRF 

m 

$ 

(rads  where 

0  <  ®  <  2*) 
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GYRFD 

OAQ  * 

F { Q * }  (°/ sec) 

Filtered  value  of  Q 

0AQ1  - 

Q'  (°/sec) 

pitch  rate  corrupted  by 

measurement  errors 

OAR 

F{  R* }  (°/sec) 

Filtered  value  of  R' 

PARI  - 

R'  (°/sec) 

yaw  rate  corrupted  by 

measurement  errors 

determined  by  the  call 

to  GYERR. 

This  subroutine  is  called  from  FLTDYN  every  INTRTE  seconds. 


Subroutine  KINEM  (VELP) 
Program  GDDRONE 


This  kinematics  subroutine  takes  linear  velocities  in  the  vehicle  frame 
and  transforms  them  to  linear  velocities  in  the  inertial  frame  using  the  Euler 
angle  defined  matrix.  These  velocities  are  integrated  to  determine  the  up¬ 
dated  vehicle  position  in  space.  The  inputs  and  outputs  are: 


INPUTS 

ARGUMENT  LIST 

VELP(3 ) 

VELP(l) 

VELP(2 ) 

VELP(  1> ) 

Vx  velocities  in 

Vy  vehicle  reference 

Vz  f  rame 

OUTPUTS 

LABELLED  COMMON 

H0J2 

P0SITN(3)  - 

3-component  vector  describing  position 

in  iner< ial  space 
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VELCTY(3) 


3-coaponent  vector  describing  velocity 
in  interial  space 


DATRNS  VN  - 

GAMPN  - 


GAMYN  - 

This  subroutine  is  called  from  FLTDYN  every  INRTE  seconds. 


Subroutine  WIND  (A,B,C,D,E,F) 
Program  GDDRONE 


This  subroutine  provides  the  wind  linear  and  angular  velocities  in  the 
vehicle  frame  of  reference.  The  arguments  are: 


V 


A 

m 

vwx  (ft/sec) 

linear  and 

B 

m 

Vwy  (ft/sec) 

angular 

C 

m 

Vwz  (ft/sec) 

> 

velocities 

D 

m 

Pw  (rads/sec) 

associated 

E 

m 

(rads/sec) 

with  the  wind 

F 

CB 

Ry,  (rads/sec) 

✓ 

field 

It  is  first  called  from  FLTDYN  for  initialization  and  then  it  is  called 
from  EQAMQ  every  1NTRTE  seconds. 


f 


l  VELCTY  (I) 


J-l 


tan 


VELCTY(  3) 


VELCTY2 (1)  +  VELCTY2 (2) 


tan 


■i  / VELCTY  (  2)\ 
\ VELCTY  (1 )/ 
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Subroutine  RCWINU  (A,8,C,D,E,F) 

Program  GDDRONE 

This  subroutine  transforms  wind  field  velocities  from  inertially  refer 
enced  coordinates  to  vehicle  referenced  coordinates.  The  arguments  are 

A  »  X  inertially  referenced  wind  component 
B  =  Y  Inertially  referenced  wind  component 
C  =  Z  inertially  referenced  wind  component 
D  =  X  body  referenced  wind  component 
E  =  Y  body  referenced  wind  component 
F  =  Z  body  referenced  wind  component 

It  is  called  from  subroutine  WIND  every  INTRTE  seconds. 


Subroutine  RCVaNG  (A,B,C,D,E,F,G,H) 
Program  GDDRONE 


This  subroutine  uses  inertially  referenced  coordinates  to  determine  azi 
muth  and  elevation  planar  angles  relative  to  the  vehicle.  The  arguments  are: 


A  ** 
B  - 

C  - 
D  ** 

E  - 
F  - 
G  » 
H  - 


X 

Y 

Z 


vector  coordinates  to  a 
point  in  space  in  iner¬ 
tial  reference  frame. 


dummy 

dummy 

dummy 

planar  azimuth  angle  relative  to  vehicle  boresight 
planar  elevation  angle  relative  to  vehicle  boresight 


It  is  called  from  subroutine  RCWIND  every  INTRTE  seconds. 
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A  -  Number  of  bits  used  to  represent  variables 
B  “  Scaling  parameters 
C  -  The  truncated  quantity 

D  -  Integer  representation  of  the  truncated  quantity 
This  subroutine  is  called  from  CRSGID  every  DRATE  seconds. 


This  subroutine  is  a  two-dimensional  array  lookup  routine.  The  arguments 

are: 

A  *  Table  1  (D  elements) 

B  »  Table  2  (E  elements) 

C  -  Table  3  (D  *  E  elements) 

F  -  Current  A  value 
G  ■  Current  B  value 
H  ■  Interpolated  value 

This  subroutine  is  called  form  AERDYN  every  INTKTE  seconds. 
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Subroutine  BLIHIT  (A,B,C) 
Program  GDDRONE 


This  subroutine  limits  a  variable  between  upper  and  lower  limits.  The 
arguments  are: 


A  «  Variable  to  be  limited 
B  -  Lower  limit 
C  »  Upper  limit 

This  subroutine  is  called  from  both  EQAMO  and  K1NEM  every  INTRTE  seconds. 


Subroutine  INTGRT  (A,B,C,D,E,F) 
Program  GDDRONE 


This  subroutine  is  the  Integration  procedure  used  in  the  program: 

A  -  A  +  D(2B  -  C) 

where  A  is  the  variable  being  integrated 

D  is  the  integration  interval 

B  is  the  value  of  A  at  the  previous  interval  (n  -  i) 
C  is  the  value  of  A  two  Intervals  prior  (n  -  2) 

E  -  1 
F  -  1 

This  subroutine  is  called  every  INTRTE  seconds  by  both  EQAMO  and  KXNEM. 
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Subroutine  INVERT  (A,B,C) 
Program  GDDRONE 


This  subroutine  transforms  a  vector  from  body  to  inertial  coordinates. 
The  arguments  are: 

A  ■  Boresight  vector  in  inertial  coordinates 
B  ■  Vector  in  body  coordinate  frame 
C  -  Vector  in  space  coordinate  frame 

This  subroutine  is  called  every  INTRTE  seconds  from  the  subroutine  KINEM. 


Subroutine  GYERR  (A.B) 
Program  GDGUID 


This  subroutine  computes  a  bias  corrupted  value  of  the  body  rates  to 
simulate  gyro  error.  The  arguments  are: 

A  »  Uncorrupted  value  of  the  body  rate  as 
input  and  corrupted  value  upon  return 
B  -  Index  which  refers  to  particular  value 
of  error  which  are  stored  in  a  common 
array. 

1  “  yaw 

2  -  dummy 

3  "  pitch 

The  values  for  yaw  and  pitch  bias  are  entered  from  the  terminal  at  run  time. 
This  subroutine  is  called  from  EC^MO  every  INTRTE  seconds. 


J 
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